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I Abstract 

(D 
rG^ The matter spectrum of the MSSM, including three right-handed 

neutrino supermultiplets and one pair of Higgs-Higgs conjugate super- 
^^ fields, can be obtained by compactifying the Eg x Eg heterotic string 

^vj and M-theory on Calabi-Yau manifolds with specific 5'[/(4) vector 

On bundles. These theories have the standard model gauge group aug- 

J^ mented by an additional gauged U{1)b-l- Their minimal content 

requires that the B-L gauge symmetry be spontaneously broken by 

a vacuum expectation value of at least one right-handed sneutrino. 

O III previous papers, we presented the results of a quasi-analytic renor- 

^^ malization group analysis showing that B-L gauge symmetry is indeed 

^ radiatively broken with an appropriate i?-L/electroweak hierarchy. In 

k> this paper, we extend these results by 1) enlarging the initial pa- 

^ rameter space and 2) explicitly calculating the renormalization group 

equations numerically. The regions of the initial parameter space lead- 
ing to realistic vacua are presented and the B-L/electroweak hierar- 
chy computed over these regimes. At representative points, the mass 
spectrum for all sparticles and Higgs fields is calculated and shown 
to be consistent with present experimental bounds. Some fundamen- 
tal phenomenological signatures of a non-zero right-handed sneutrino 
expectation value are discussed, particularly the cosmology and pro- 
ton lifetime arising from induced lepton and baryon number violating 
interactions. 
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1 Introduction 

Es X Eg heterotic strings [Ij and heterotic M-theory ^Sj-lZj offer perhaps the 
simplest approach for deriving reahstic particle physics from superstrings. 
Their compactification on Calabi-Yau threefolds with slope-stable holomor- 
phic vector bundles leads to A^ = 1 supersymmetric effective theories in 
four-dimensions. Such vacua include complete intersection and elliptically 
fibered Calabi-Yau spaces admitting vector bundles constructed using mon- 
ads [8]- [12], spectral covers [13]- [16] and extension of lower rank bundles 
[T71 HH]. The formalism for computing the low energy spectrum in each 
case has been developed, and presented in [IHl [20], [211 ES] and [231 El] re- 
spectively. Cohomological methods have been used to calculate the texture 
of Yukawa couplings and other parameters in these contexts [25]- [27]. Fi- 
nally, non-perturbative string instanton contributions to the superpotential 
are computed in [2H]-[3I] and used to discuss moduli stability, supersymme- 
try breaking and the cosmological constant [32] • These methods underlie the 
theory of "brane universes" [HI |33] . 

For certain choices of Calabi-Yau manifolds with non-trivial homotopy, 
along with vector bundles with the appropriate structure group and topol- 
ogy, the derived low energy theory can be phenomenologically viable |34j - 
|37j . Specifically, for Calabi-Yau manifolds with Z3 x Z3 homotopy and 
a vector bundle with S't/(4) structure group, it has been shown [3H1 [39] 
that the low-energy theory can have exactly the matter and Higgs spec- 
trum of the MSSM, including three right-handed neutrino supermultiplets, 
one per family, along with a relatively small number of uncharged geomet- 
ric and vector bundle moduli. The gauge group of this effective theory is 
SU{3)c X SU{2)l X U{1)y x U{1)b-l, that is, the standard model gauge 
group times an additional B-L Abelian gauge symmetry. We will refer to 
this as the B-L MSSM theory. 

The existence of the extra U{1)b-l gauge factor, far from being being ex- 
traneous or problematical, is precisely what is required to make a heterotic 
vacuum with SU{4) structure group phenomenologically viable. The reason 
is the following. As is well-known, four-dimensional A^ = 1 supersymmet- 
ric theories generically contain two lepton number violating and one baryon 
number violating dimension four operators in the superpotential. The former, 
if too large, can create serious cosmological difficulties, such as in baryogen- 
esis and primordial nucleosynthesis [10]- [13], as well as coming into conflict 
with direct measurements of lepton violating decays [B] . The latter can pro- 



duce extremely rapid proton decay, far in excess of the observed bound on 
its lifetime [^ H5] . To avoid these problems, it is traditional in low-energy 
N = 1 supersymmetric theories to impose a discrete "matter parity", the 
supersymetric version of "R-parity" |15]-[5T]. This Z2 finite symmetry dis- 
allows these dimension four operators from appearing in the superpotential, 
thus solving all the above problems. The B-L MSSM theory naturally con- 
tains matter parity as a Z2 subgoup of U{1)b~l- As long as this subgroup 
is unbroken, or weakly broken, the theory will be phenomenologically and 
cosmologically viable. Importantly, however, since a gauged B-L Abelian 
symmetry is not observed at low energy, it is essential that U{1)b-l be 
spontaneously broken above the electroweak scale. 

In pre-superstring supersymmetric grand unified theories (GUTs), this 
phenomenon was discussed within the context of the unification group 5*0(10) 
[52]-[57j. It was noted that certain SU{3)c x SU{2)l x U{1)y singlets carry- 
ing non-zero 3(i? — L) even charge would, if they obtained a non-vanishing 
vacuum expectation value (VEV), spontaneously break U{1)b-l down to a 
residual Z2 symmetry, that is, matter parity. The superpotential is then con- 
structed so that their VEVs are very large, usually near the unification scale. 
Thus, these multiplets are heavy and disappear from low-energy physics, 
leaving behind unbroken matter parity. However, as first pointed out in [58] , 
this compelling mechanism cannot occur in smooth heterotic compactifica- 
tions. The reason is straightforward. The requisite 3{B — L) even charged 
singlets can only arise in the decomposition of 5*0(10) representations with 
dimension 126 and higher. Although unconstrained GUTs can always add 
such multiplets, they can never appear as zero-modes of the Dirac opera- 
tor in smooth 5f/(4) heterotic compactifications. All such zero-modes must 
arise from the decomposition of the 248 representation of Eg under 5f/(4), 
whose largest dimensional 5*0(10) representation is the 45. It follows that 
all SU{3)c X SU{2)l x U{1)y x U{1)b-l muhiplets that appear after Wil- 
son line breaking will have charges 3{B — L) = ±1, ±3, 0. Hence, the above 
mechanism to obtain matter parity never occurs. 

What, then, can one do in smooth heterotic 5*f/(4) compactifications? 
Note that the only 5*f/(3)c x SU{2)l x f/(l)y singlets carrying non-zero 
U{1)b-l charge are the right-handed sneutrinos, each of which has 3{B—L) = 
3. Were at least one sneutrino to develop a non-zero VEV, it would sponta- 
neously break the U{1)b~l symmetry. However, since its B-L charge is odd, 
the Z2 matter parity will also be broken. It follows that a viable theory can 
only emerge if at least one right-handed sneutrino develops a non-vanishing 



VEV that is 1) larger than the electroweak scale but 2) small enough that 
the broken matter parity can sufficiently prevent large lepton and baryon 
number violation. That is, one should have a reasonable i?-L/electroweak 
hierarchy. It was shown in [581 EH] that this radiative hierarchy can, indeed, 
occur. Specifically, using a quasi-analytic solution to the renormalization 
group equations (RGEs) it was found that a non-zero VEV can develop in 
any right-handed sneutrino followed, at a sufficiently lower scale, by the stan- 
dard radiative breaking of electroweak symmetry via Higgs VEVs. That is, 
the B-L MSSM theory can be a viable theory of nature with interesting cos- 
mological consequences [20] . However, the analysis in [SHI EH] was restricted 
in two important ways. First, in order to attain a quasi-analytic solution it 
was necessary to choose realistic, but constrained, initial data for the soft 
supersymmetry breaking parameters. Second, the quasi-analytic solution re- 
quired that relatively "small" terms were dropped in the calculation. Hence, 
for example, the superparticle mass spectra presented in [581 EH] were only 
approximate, as were the calculations of the 5-L/electroweak hierarchy. 

In this paper, we present a completely numerical calculation of the RGEs 
in the B-L MSSM theory. Hypercharge and B-L kinetic mixing does occur 
in this theory and, generically, should be included in the analysis. However, 
we show in Appendix C that for the regions of parameter space of emphasis 
here, this mixing is small and, to leading order, can be ignored. This allows 
a great simplification in a theory which already has a large number of input 
parameters. With this caveat, all RGEs are solved without approximation. 
Furthermore, one can now explore the complete initial soft supersymmetry 
breaking parameter space and find the sub-regions that give realistic particle 
physics and cosmology. Specifically, we do the following. In Section 2, the ex- 
act MSSM spectrum and the associated SU{3)c x SU{2)l x U{1)y x U{1)b-l 
quantum numbers are presented, along with the superpotential, D-terms and 
soft supersymmetry breaking quadratic and cubic terms. Section 3 is devoted 
to extending the ideas in [58l [59] relevant to ensuring the spontaneous break- 
ing of U{1)b-l through right-handed sneutrino VEVs, as well as specifying 
some physically less interesting initial parameters. The number of initial 
parameters is reduced to four, related to the squarks, right-handed sneutri- 
nos, the fi parameter and tan/3. Three phenomenological constraints are 
then presented. The first is two inequalities that ensure radiative breaking of 
electroweak symmetry through up- and down-Higgs VEVs. Second, we give 
the constraints required to make the i?-L/electroweak vacua local minima 
of the potential energy. Third, the lower bounds on the masses of all super- 



partners, as well as the Higgs fields, are presented. The calculations in this 
paper will satisfy all three constraints. Our main numerical results are pre- 
sented in Section 4. This section is broken into three subsections and a brief 
summary. The three subsections reflect the fact that all B-L MSSM vacua 
can be catagorized by the sign of the left-handed squark and right-handed 
slepton squared masses; that is, 1) all rn? > 0, 2) all masses positive except 
rriQ^ < and 3) all masses positive except m^. < 0. In each case, we find the 
complete region of parameter space for which one obtains a realistic theory, 
compute the S-L/electroweak hierarchy over each acceptable region and, at 
some representative points, explicitly compute the sparticle and Higgs mass 
spectrum. We verify that phenomenological mass constraints are indeed sat- 
isfied at these points. A discussion of the formalism used in this paper to 
calculate fermion and scalar masses is presented in Appendix A. The present 
experimental constraints on the Higgs masses are reviewed in Appendix B. 

Our results are predictive, since many low-energy phenomena arise from 
the radiative breaking of a right-handed sneutrino. Perhaps the most strik- 
ing aspect of this is that the non-vanishing sneutrino VEV "grows back" 
the previously disallowed lepton number violating dimension four terms in 
the superpotential, each with an explicitly calculable coefficient. Following 
[GTI [62] . we confront our results with various cosmological constraints, such 
as baryon asymmetry and primordial nucleosynthesis. We find that they are 
all easily satisfied in the B-L MSSM theory. Furthermore, we show that our 
theory is consistent with gravitino dark matter and rapidly decaying stan- 
dard model sparticles. Another important aspect of breaking B-L symmetry 
with a right-handed sneutrino is that the previously disallowed baryon vio- 
lating dimension four operator does not grow back from the dimension four 
superpotential. It can only reappear from higher dimensional operators with 
calculable, and naturally suppressed, coefficients. Putting in our calculated 
results, we find that proton decay through dimension four operators can be 
sufficiently suppressed to satisfy all bounds on the proton lifetime. These 
lepton and baryon number violating results are presented in Section 5. 

Finally, we want to point out that the B-L MSSM theory, with sponta- 
neous breaking of U{1)b-l through right-handed sneutrinos, was presented 
from a "bottom up" point of view in [63]-[65j. These authors discussed 
various phenomenological predictions and applied similar ideas to other low- 
energy theories [66] . 



2 The A^ = 1 Supersyminetric Theory 

We will consider an A^ = 1 supersymmetric theory with gauge group 

G = SUi3)c X SUi2)L X f/(l)y X f/(l)B-L (1) 

and the associated vector superfields. The gauge parameters are denoted 
by Qs, 92, Qy and Qb-l respectively. The matter spectrum consists of three 
families of quark and lepton chiral superfields, each family with a right-handed 
neutrino. They transform under the gauge group in the standard manner as 

g, = (3,2,1/3,1/3), M, = (3,1,-4/3,-1/3), rf, = (3, 1, 2/3, -1/3) (2) 

for the left and right-handed quarks and 

L, = (1,2, -1,-1), z/, = (l,l,0,l), e, = (l,l,2,l) (3) 

for the left and right-handed leptons, where i = 1, 2, 3. In addition, the spec- 
trum has one pair of Higgs-Higgs conjugate chiral superfields transforming 
as 

iJ= (1,2,1,0), i? = (l,2,-l,0). (4) 

When necessary, the left-handed SU{2)l doublets will be written as 

Q, = (f/„A), L, = {N„E,), H={H+,H''), H={H\H-). (5) 

There are no other fields in the spectrum. 

The supersymmetric potential energy is given by the usual sum over the 
modulus squared of the F and D-terms. In principle, the F-terms are de- 
termined from the most general superpotential invariant under the gauge 
group, 

3 

W = fj,HH + ^ {\u,ijQiHuj + \d,ijQiHdj + X^^ijLiHuj + Xe^ijLiHcj) (6) 

Note that the quadradic mixing term of the form LiH, as well as the dan- 
gerous lepton and baryon number violating interactions 



which generically would lead, for example, to rapid nucleon decay, are disal- 
lowed by the U{1)b-l gauge symmetry. To simplify the upcoming calcula- 
tions, we will assume that we are in a mass-diagonal basis where 

^u,ij = ^d,ij = ^u,ij = ^e,ij = 0, i ^ j. (8) 

Note that once these off-diagonal couplings vanish just below the compacti- 
fication scale, they will do so at all lower energy-momenta. We will denote 
the diagonal Yukawa couplings by Xa = Xi, i = 1,2,3. Next, observe that 
a constant, field-independent fi parameter cannot arise in a supersymmet- 
ric string vacuum since the Higgs fields are zero modes. However, the HH 
bilinear can have higher-dimensional couplings to moduli through both holo- 
morphic and non-holomorphic interactions in the superpotential and Kahler 
potential respectively. When moduli acquire VEVs due to non-perturbative 
effects, these can induce non-vanishing supersymmetric contributions to fi. 
A non-zero yU can also be generated by gaugino condensation in the hidden 
sector. Why this induced /i-term should be small enough to be consistent 
with electroweak symmetry breaking is a difficult, model dependent problem. 
In this paper, we will not discuss this "yU-problem" . Instead, we will consider 
the /i parameter as an input to our analysis and consider a range of possible 
values. 

The SU{3)c and SU{2)l -D-terms are of the standard form. We present 
the U{1)y and U{1)b-l -D-terms, 

Dy = ^y + 9y<P\{Y/2)^^<Pb (9) 

and 

Db-L = iB-L + gB~L4>\ (Yb-l)ab 0B (10) 

where the index A runs over all scalar fields (pAi to set the notation for the 
hypercharge and B-L charge generators and to remind the reader that each 
of these D-terms potentially has a Fayet-Iliopoulos (FI) additive constant. 
However, as with the /i parameter, constant field-independent FI terms can- 
not occur in string vacua since the low energy fields are zero modes. Field- 
dependent FI terms can occur in some contexts, see for example [HZj- How- 
ever, since both the hypercharge and B-L gauge symmetries are anomaly 
free, such field-dependent FI terms are not generated in the supersymmetric 
effective theory. We include them in (|9|,(10) since they can, in principle. 



arise at a lower scale from radiative corrections once supersymmetry is softly 
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broken [BH]- Be that as it may, if calculations are done in the D -eliminated 
formalism, which we use in this paper, these FI parameters can be consis- 
tently absorbed into the definition of the soft scalar masses and their beta 
functions. Hence, we will no longer consider them. 

In addition to the supersymmetric potential, the Lagrangian density also 
contains explicit "soft" supersymmetry violating terms. These arise from the 
spontaneous breaking of supersymmetry in a hidden sector that has been inte- 
grated out of the theory. This breaking can occur in either F-terms, D-terms 
or both in the hidden sector. In this paper, for simplicity, we will restrict 
our discussion to soft supersymmetry breaking terms arising exclusively from 
F-terms. The form of these terms is well-known and, in the present context, 
given by [69]-[73] 

l^soft = l"2s + Vss + V2f, (11) 

where V2s are scalar mass terms 

3 

TA \ ^/ 2 |,n |2 I 2 I |2 I 2 I J |2 I 2 I r |2 I 2 i |2 

V2s = 2^{mQ^\Qi\ +m^JUi\ +m^^\di\ +m^JLi| +m^Ji/i| 

+ml\ei\^) + mjj\H\'^ + m}j\H\^ - {BHH + he), (12) 

Vss are scalar cubic couplings 

3 

Vss = J2(^^iQiHui + Ad^QiHdi + A^XiHui + Ae^Ha + he) (13) 
and V2/ contains the gaugino mass terms 

V2f = -M3A3A3 + -M2A2A2 + -MyAyAy + -Mb-L^B-lXb-L + hc. (14) 



As above, to simplify the calculation we assume the parameters in (12) and 



(13) are flavor-diagonal. This is consistent since once the off-diagonal param- 
eters vanish just below the compactification scale, they will do so at all lower 
energy-momenta. Finally, note that lepton and baryon violating scalar cubic 
terms of the form ([T]) are disallowed in V^s by the U{1)b-l gauge symmetry. 

3 Initial Parameter Space 

The four- dimensional effective theory described in the previous section arises 
at an initial energy-momentum just below the compactification scale given by 



the inverse Calabi-Yau radius. In order to carry out a detailed renormaliza- 
tion group analysis, we must specify this initial energy-momentum precisely. 
We will do this as follows. 

3.1 Gauge Coupling Parameters 

It is well known that precision measurements [71]- [TB] carried out at the 

electroweak scale indicate that the SU{3)c x SU{2)lX f/(l)y gauge couplings, 

g3,g2 and gi = y fs'y respectively, unify to 

g{0) ~ .726 (15) 

at scale 

M„ ~ 3 X lO^^GeV . (16) 

For simplicty, so that we can ignore a discussion of threshold effects, we will 
assume that the initial energy momentum for our effective theory is precisely 
the unification scale M„. In addition, since the S'[/(4) vector bundle breaks 

Es to 50(10), we will take the U{1)b-l gauge coupling g^ = J^gB-i to 
unify with the three other couplings at M„. 

Having fixed the initial energy- momentum as M„, one must now specify the 
initial values of all parameters in the effective theory at this scale. In princi- 
ple, string theory would predict these parameters as functions of the moduli 
VEVs. In this paper, however, we will be content with simply choosing the 
initial parameters subject to the dictates of simplicity, the "universality" of 
some parameters observed in minimal supergravity and simple string com- 
pactifiacations [731 177] and the necessity to break U{1)b-l through a VEV of 
at least one right-handed sneutrino. Having chosen all the initial parameters, 
their values at any lower scale, specified by 

are determined by the associated renormalization group equations (RGEs). 
These are discussed in detail in several reviews, see, for example [68]- [82], and 
were generalized to include the U{1)b-l symmetry in our previous papers 
[55| I5U] . In this paper, all calculations will be carried out at the one- loop 
level. 
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The initial unified gauge coupling is given in (15). We now turn to spec- 
ifying the initial values for all other parameters in our effective low-energy 
theory. We begin with the dimensionful parameters. 

3.2 Gaugino Mass Parameters 

Consider the soft supersymmetry breaking gaugino mass parameters that 



appear in V2/ in (14). Following standard notation, we henceforth denote 
My = Ml and Mb-l = M4. We now make the assumption that at the 
compactification scale the gaugino masses unify, that is, 

|Mi(0)| = |M2(0)| = |M3(0)| = |M4(0)| . (18) 

Such universal gaugino masses naturally occur in minimal supergravity |70j - 



and simple string theories jTTl [72] . Here, we choose (18) for reasons of 
simplicity. 

3.3 Higgs, Squark and Slepton Masses 

The RGEs for the soft supersymmetry breaking Higgs, squark and slepton 
masses all contain a term proportional to gfS where 

Y 

S = Tr(— m^) (19) 



2 2 I Y^/ 2 o 2 I 2 2,2 



It greatly simplifies the boundary condtions of these RGEs to choose the 
initial soft breaking masses so that 5(0) =0. A natural way to achieve this 
is to impose a separate unification of the Higgs masses, squark masses and 
the left doublet/down right singlet slepton masses. That is, we henceforth 
choose 

mni^f = rn^i^f, mQ^^f = m„,(0)' = m^M^ (20) 

and 

mLM' = meM (21) 

for all 2, j, k = 1,2, 3. In addition to the hypercharge induced gfS term, the 
gauged U{1)b-l symmetry of our effective theory introduces a new term into 



the RGEs for the squarks and slepton soft supersyminetry breaking masses. 
This term is of the form g^S' where 



S' = Tr{YB-Lin^)=S'o + S[ 



(22) 



and 



S'o = E(2"^Q. - < - < - < + O' '^1 = E(-< + ^'^>)- (23) 



i=l 



j=l 



It follows from (20) and (21) that Sq{0) = 0. Note, however, that unlike 
S and S'q, the S[ term depends on the soft supersymmetry breaking right- 
handed sneutrino masses. We choose the initial values of these parameters 
not to be degenerate with the other slepton masses, that is. 



mL,(0)2 = me^(0)V^..(0)' 
for all i,j, k = 1,2, 3. It follows that 



(24) 



5((0) = ^(-m^.(0)2 + m,(0)^)^0. 



(25) 



j=i 



This asymmetry is an important ingredient in generating radiative breaking 
of the U{1)b-l symmetry. We point out that soft scalar masses need not be 
"universal" in string theories, since they are not generically "minimal" . 

3.4 The A and B Parameters 

Now consider the soft supersymmetry breaking up/down Ai and B parame- 



ters in equations (13) and (12) respectively. As already stated, we take the 



Ai coefficients to be flavor diagonal. In addition, it is conventional [73j to let 



A, 



Xu.Au^, Ad^ — \d,Ad^, A 



K^Ay^, Ae^ 



KAe 



(26) 



for i = 1,2,3, where Aj are the Yukawa couplings and the dimensionful Ai 
parameters are chosen to be of order the supersymmetry breaking scale. This 
is not a requirement in the "non-minimal" string vacua that we are discussing. 



Be that as it may, for simplicity of presentation we will assume (26) for the 



remainder of this paper. The input Yukawa parameters will be discussed 
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below. In this paper, we will, for simplicity, assume the Ai parameters unify 
at the scale M„. That is, 

i„,(0) = i,^(0) = i.,(0) = ie,(0) (27) 

for all i,j, k,l = 1, 2, 3. 

The initial value of the soft breaking B parameter, -B(O), is taken to be 
arbitrary. However, in our analysis B will be treated differently than the other 
dimensionful parameters. As will be shown below, rather than choosing the 
value of the B parameter, we will instead input tanP and the supersymmetry 
breaking scale. This will dynamically fix the value of B for any given set of 
initial conditions. 



3.5 The fi Parameter 

The supersymetric fi parameter has a fundamentally different origin than 
the soft supersymmetry breaking dimensionful couplings discussed above. In 
this paper, we will simply allow its initial value /i(0) to be arbitrary. As in 
conventional radiative breaking scenarios, to be compatible with electroweak 
symmetry breaking we expect it to be of O{100)GeV. However, we make no 
attempt to solve this "/x-problem" . Having discussed the initial values for 
the dimensionful parameters, we now consider the dimensionless parameters 
in our effective theory. 

3.6 Tan/3 and the Yukawa CoupUngs 

As with any MSSM-like model, our low energy theory requires two Higgs chi- 
ral supermultiplets, H and H, whose VEVs (H) and (H) break electroweak 
symmetry and give mass to the W^ and Z vector bosons. The experimentally 
measured vector boson masses put a constraint on these VEVs. In terms of 
the Z mass, this is 

{H^ + {Hf = -^, ^ &eVr . (28) 

Hence, giving one Higgs VEV completely determines the other. It is conven- 
tional to re-express the remaining Higgs VEV in terms of the ratio 

tan/3 = H|. (29) 
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If the value of tan/3 is given, one can easily find both Higgs VEVs using (28) 



and (29). The result is 



W ^ (^Gen^iS£= . iH) . i'-^GeV)^=l=. (30) 

V2 Vl + tan2/3 V2 Vl + tan^/J 

In this paper, we will take tan/3 as an input parameter 

In addition to the vector bosons, the Higgs VEVs {H) and {H) give mass 
to the up and down quarks/leptons respectively. As with the SU{3)c x 
SU{2)l X f/(l)y gauge couplings, the Yukawa couplings are highly con- 
strained by experiment. Given a value of tan/3 and, hence, {H) and (H), 
the known masses of the quarks/leptons completely determine the Yukawa 
couplings at the electroweak scale. However, unlike the gauge couplings, the 
Yukawa coupling do not unify at M„. Rather, when run up to the unification 
scale using their RGEs, the initial values of the Yukawa couplings are a set 
of tan/3 dependent numbers with no particular relationship. Therefore, in 
this paper, rather than specifying the initial Yukawa couplings at scale M^, 
we will instead input a value of tan/3 and use the associated Higgs VEVs 
and the measured quark/lepton masses to calculate all Yukawa parameters 
at the electroweak scale. These will then be run back to the unification scale 
and stored in our program. When required, the initial Yukawa parameters 
can then be input into any other ROE and scaled down along with the other 
relevant parameters. 



It is important to note from (30) that as tan/3 is decreased, the up Higgs 
VEV (H) must get smaller. This then necessitates taking larger values for 
the up Yukawa couplings to be consistent with the measured masses. For (H) 
sufficiently small, the top quark Yukawa coupling will become much larger 
than unity and the theory becomes non-perturbative. This puts a bound 
on how small (H) can be and, hence, a lower bound on tan/3. Similarly, 
increasing tan /3 requires the down Higgs VEV {H) to decrease. For {H) suf- 
ficiently small, the bottom quark Yukawa coupling will become much larger 
than unity and the theory non-perturbative. This puts a bound on how small 
{H) can be and, hence, an upper bound on tan/3. These bounds on tan/3 
are typically estimated [75], |H5] to be 

4 < tan/3 < 50. (31) 

When inputting tan/3 in this paper, we will always restrict it to be within 
these bounds. 
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3.7 Parameterizing the Initial Conditions 

Recall that, with the exception of the fi parameter, all of the dimensional 
coefficients discussed above occur in soft supersymmetry breaking interac- 
tions. If we denote by A^ a mass characterizing the scale of supersymmetry 
breaking, then each of the above coefficients can be written in the form 

Ciit)M , (32) 

where Cj(t) is dimensionless. This parameterization emphasizes that the soft 
dimensionful coefficients share a common supersymmetry breaking scale. The 
initial coefficients, Cj(0), are arbitrary. However, naturalness would dictate 
that they not to be too much larger, or smaller, than unity. The exception to 
this is the parameter fi. This arises in the supersymmetric quadratic Higgs 
term and is, a priori, unrelated to the scale Ai. However, it can always 



be written in the form (32). In this case, however, one does not expect the 
associated coefficient to be of order unity. Be that as it may, the "/^-problem" 
specifies that appropriate radiative electroweak breaking will require fi to be 
close to the scale JH. 

Specifically, this parameterization of the dimensionful parameters allows 
us to write the initial value for the gaugino masses as 

|Mi(0)| = |M2(0)| = |M3(0)| = |M4(0)| = cm{0)M , (33) 

as well as 

mniO) = ruHiO) = CHiO)M (34) 

for the initial Higgs parameters. Similarly, the initial squark and doublet/ 
down singlet slepton masses are 

mQXO) = m„^.(0) = rridM = Cg{0)M (35) 

and 

rriLXO) = m.^iO) = CeiO)M (36) 

respectively. However, for the reasons discussed below, we will allow the 
initial right-handed sneutrino masses to have the texture 

m,, (0) = m,, (0) = c,,_, {0)M , m,.^ (0) = 0,3 {0)M . (37) 

Finally, we write 

i.,(0) = i,^.(0) = i.,(0) = 4(0) = c^iO)M (38) 
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and 

/i(0) = c^iO)M , 5(0) = cl{0)M^ (39) 

for the initial dimension-one A, fi parameters and the dimension-two B pa- 
rameter respectively. That is, there is a total of nine dimensionless Cj(0) 
parameters arising from the dimensionful parameters in our effective theory. 
However, this number can be reduced as follows. 

First, note that all mass parameters scale with the same factor Ai. Hence, 
one can always redefine Ai so as to absorb one of these coefficients. Without 
loss of generality, we can choose this to be the Higgs parameter. That is, set 

cniO) = 1 . (40) 

Second, in minimal supergravity and simple superstring vacua, the unified 
initial A and gaugino mass parameters are numbers of order unity times the 
supersymmetry breaking scale Ai. We will assume this in our calculation as 
well. For simplicity, choose 

Ca(0) = 1 • (41) 

The initial value for cm is more subtle to determine. We have done an exten- 
sive numerical analysis of phenomenologically acceptable initial conditions al- 
lowing cm(0) to vary freely. The result is a bound given by 0.1 < cm(0) < 1.2. 
In this paper, for simplicity of presentation, we fix this initial parameter to 
a value in the middle of this range given by 

cm(0) = 0.6 . (42) 

Finally, we will also specify the coefficients Ce(0) and 0,^^2(0) as follows. 

In a previous paper [59], we presented a quasi- analytic solution to the 
RGBs in the B-L MSSM theory subject to certain initial conditions on the 
parameters. To obtain an analytic solution, the initial parameters chosen 
were considerably more constrained than they are in this paper. Be that as 
it may, the generalized parameter space discussed here contains these initial 
conditions as a small subset. Specifically, we showed that at the B-L scale 
Mb-l — lO^GeT^ the right-handed down slepton and right-handed sneutrino 
soft mass parameters are given by 

me, {tB-L? = m,, {Of - (3.35 x 10-^)S[ (0) , (43) 

m.^tB-Lf = m,,{Oy - (3.35 x 10-^)S[{0) (44) 
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for 2 = 1,2,3 where, using (25), (36) and (37), one can write 

with 

For specificity, let us choose 

(3.35 X 10"2)(1 + 2^2 - 3^2) = 5 . 
Then one obtains the simple result that 



my.j^{t 



B-L 



-4m,3(0)^ 



(45) 
(46) 

(47) 
(48) 



leading to a non-zero VEV in the u^ direction. In this way, we guarantee 



radiative U{1)b-l breaking in the theory. Similarly, using (47) we find that 

m,,,, (tB-LY = (C^ - 5)m,3 (0)2 (49) 

and 

me,(tB_i)' = (A2-5)m,3(0)2 (50) 

for i = 1,2, 3. The simplest vacuum structure occurs when all meXtB-hY are 
positive. For this to be the case, the coefficient A must satisfy A^ — 5 > 0. 
Again, for specificity we will choose 

A = VQ , (51) 

which yields the simple result that 

meXtB-bf = my^{Qif 

for 2 = 1, 2, 3. Putting A = -\/6 into expression (47) gives 

C = 9.12 . 

It then follows from (49) that both mu^ 2(^-8-^)^ ^^^ positive and given by 

(54) 



(52) 



(53) 



m,y,^Xt 



B-L 



78.2 m^3(0)^ . 



We conclude that the choice of the A and C parameters given in (51) and 



(53 ) respectively leads to a vacuum that has positive soft squared masses and. 
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hence, vanishing VEVs for all sleptons with the exception of the third family 
right-handed sneutrino. This acquires a non-zero VEV which radiatively 
breaks U{1)b-l symmetry. It is clear that these choices for A and C are far 
from unique, and that a wide range of values would still lead to a vacuum 
with appropriate U{1)b-l symmetry breaking. Be that as it may, we find 



it convenient to continue to use (51) and (53) in the present paper. It then 



follows from (46) that we will choose 



Ce{0) = V6 c,, (0) , c,,,, (0) = 9. 12 c,3 (0) 



(55) 



The constraints given in (^40h, (41), (42) and (55) reduce the number of 



free parameters down to six- four Cj parameters as well as Ai and tan/3. 
There are, however, important phenomenological constraints on these pa- 
rameters, to which we now turn. 



3.8 Phenomenological Constraints 

It is well-known [73] that for an MSSM-like theory with two Higgs doublets, 
H and H, to have a stable vacuum solution that breaks electroweak sym- 
metry, the parameters of the theory have to satisfy two constraints at the 
electroweak scale Mew — lO'^GeV. These are 



i?2>(|/ip + m|)(|/ip + m|) 



(56) 



which ensures that one linear combination of H and H has a negative squared 
mass, thus enabling a non-zero Higgs VEV to form, and 



2B < 2|/ip + m^ + m^ 



(57) 



which guarantees that the quadratic part of the potential energy is positive 
along the D-flat directions and, hence, that the potential energy is bounded 
from below. Once these conditions are satisfied, the theory has a stable Higgs 
vacuum specified by the two minimization equations. Their solutions can be 
put in the form 

sin(2/3) = , , ^f , „ „ (58) 



771 



H 



+ 7njy + 2\n\ 



and 



Ml 



12 2 1 

\7n^ — 77ljj\ 

l\ - sin2(2/3) 



7Jl 



H 



771 



H 



-2|/i| 



(59) 
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with the parameters evaluated at the electroweak scale. 

In many analyses of electroweak breaking, all the soft masses and the jj, 



parameter are given as input with tan (3 and Mz generated as solutions of (58 ) 



and (59). However, as discussed above, it is convenient in this paper to take 



tan/3 as an input parameter. It follows that equation (58) should be viewed 
as yet another constraint on the soft breaking parameters. Specifically, we 
will use (58) to solve for 5 as a function of tan/3, rnjj, ni?fr and /i at the 



electroweak scale. This is possible since the RGEs for ?ti|^, m^ and /z [79] 
and, hence, the value of these parameters at the electroweak scale do not 
depend implicitly on B. Written in terms of the notation introduced in the 
previous section, it follows that 

2 sin(2f3) 



-B 



ic 



' +4 + 2|c. 



H 



■'H 



(60) 



We can then scale this parameter back up to the M^ to determine the initial 
value 5(0). 



Similarly, we can input the experimental value of Mz into (59) and use 



this to put a further constraint on the initial parameters. In terms of the 



above parameterization, ( 59 ) can be re-written as 



Mi 



c^-c 



H\ 



■^H 



r^ -9lr |2 



M' 



(61) 



Vl-sin2(2/3) 

From this equation we see that, given the initial values of q and tan/3, one 
can use the experimentally derived value for Mz to fix M. and, thus, the soft 
breaking scale. Note that for fixed values of ch, eg and tan/3, mass A^ is a 
minimum as c„ — )■ and becomes arbitrarily large as 



U- 



^H 



- r2 I 



-Cfj-C 



H 



(62) 



2V^l_sin2(2/3) 
It follows that the value of the supersymmetry breaking parameter is not 



particularly restricted by constraint (61). Be that as it may, obtaining its 
minimum value and, in particular, a large value requires fine-tuning c^ to 



zero and (62) respectively. Without fine-tuning, the typical value for M. is 



set by the Z-mass and, for the initial parameters in this paper, found to be 
of order a few hundred GeV up to order 10 TeV. 



Applying constraints (60) and (61) to fix the values of cb(0) and M. 
respectively, we are now left with four free parameters. They are 



c,(0), 



C^3(0), 



c^(0), tan/3 



(63) 
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In the remainder of this paper, we analyze the vacuum state and mass spec- 
trum of the B-L MSSM theory over this four-dimensional initial parameter 
space. This is accomplished by numerically solving all RGEs for a given 
choice of initial conditions, scaling down from M^ to Mew- Iii doing this, 
however, we will impose several important phenomenological constraints- 
rejecting the initial parameters if the results fail to satisfy these constraints 
and accepting them if they are satisfied. In this way, one can map out the 
allowed region of the four-dimensional initial parameter space. 

The phenomenological constraints we impose are the following. 

• To ensure that a stable electroweak breaking vacuum can develop at 



low energy- momenta, we impose the constraint that inequalities (56) 



and (57) be satisfied. This should be understood as a consistency check 
on our assumption, implicit in using tan /3 and the experimental value 
of Mz as input parameters, that a stable electroweak breaking vacuum 



described by (58) and (59) exists. In terms of the parameterization 



introduced in Subsection 3.7, these constraints are 



and 



respectively. 



4>(ic,r+4)(ic,r+4) (64) 



24<4 + c| + 2|c^|2 (65) 



As discussed above, condition (55) ensures that a vacuum expectation 
value develops in the third right-handed sneutrino. To guarantee that 
this is a stable local minimum, we impose the constraint that the effec- 
tive squared masses of all squarks and sleptons evaluated at the B-L 
breaking VEV (z/3), for example, 

(m|J = ml^ + \gl{v,f , (mij = ml^ - \gl{v,Y (66) 

are positive over the entire scaling range. It follows that color and 
charge symmetry are never spontaneously broken. Note that imposing 
the positivity of the effective masses does not necessarily restrict the 
soft squared masses to be positive. For example, the positivity of (^gj 
does not require that toq. be positive. On the other hand, m\^ must 
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lable L'. Experimental lower bounds on the Higgs fields and sparticles in the MSSM. The Z mass 
is for an additional U{1) gauge boson arising from spontaneously broken 50(10). See Appendix B for a 
discussion of the bounds on the Higgs scalars. 

be positive to ensure that («^|.) is. This allows us to classify the B-L 
MSSM vacua in terms of the signs of the soft squared masses at the 
electroweak scale. This will be discussed in detail later in the paper. 

• An important phenomenological constraint is that our results be con- 
sistent with the observed bounds on the masses of the Higgs fields, Hig- 
gsinos and all squarks, sleptons and gauginos. Note that such bounds 
are partially dependent on the theory in which they are analyzed. The 
determination of the exact experimental bounds within the context of 
the B-L MSSM theory is left to future work. In this paper, we use the 
fact that this theory is a minimal extension of the MSSM and, hence, 
expect most bounds to be similar to those found for the MSSM. The 
most conservative bounds for the MSSM are given in the Summary Ta- 
bles contained in the Particle Data Group review [75] and reproduced 
in Table [ij We emphasize that these serve as guidelines rather than 
strict bounds, since we are working with a model that is somewhat 
different than the MSSM. 

The theoretical calculation of the masses in the B-L MSSM involves 
considerable mixing of the fields induced by the z/3 and H,H VEVs. 
This presents a challenge in our analysis. Details of the mass matrices, 
the diagonalization process, as well as a discussion of the role of the 
spontaneously broken B-L gauge symmetry, are presented in Appendix 
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A. In this paper, we compute the mass eigenvalues for the Higgs fields 
and all sparticles and compare the results to the values in Table [l] 
We disallow all initial conditions that violate these bounds. This re- 
quires particular care for the Higgs fields, and is discussed in detail in 
Appendix B. 

4 Numerical Analysis 

We now turn to the numerical analysis of the low-energy vacua associated 



with the four initial parameters given in (63). Even though the number of 
these parameters has been reduced to four, a systematic study of this space 
is still labor intensive. Happily, there is a natural splitting into two two- 
dimensional spaces. To see this, note that one of the physical properties we 
are most concerned with is the hierarchy between the B-L and electroweak 
breaking. This hierarchy can be described in several ways [59]. Here, we will 
define the hierarchy as the ratio of the mass of the U{1)b-l gauge boson, 
given by 

MA,_, = V2gB-L{v,), (^3) = ^^ (67) 

Qb-l 

evaluated at the electroweak scale, and the Z-boson mass given in ( [59| ). 



Written in terms of the parameterization introduced in Subsection |3.7| the 
hierarchy becomes 



Ma,_, V2| 



_ ^ CggN 

The factor of M. occurs in both the numerator and the denominator and, 
hence, cancels out of this expression. Of the five parameters in (68), only Cy^, 



c^ and tan (3 have arbitrary initial conditions. Noting that all q coefficients, 
even when evaluated at the electroweak scale, are essentially of order unity, 
we see that the most infiuential factors in the size of the hierarchy are c^ and 



tan/3. This is because for fixed tan/3 one can drive the denominator in (68) 
to zero, and, hence, the hierarchy to be arbitrarily large, by fine-tuning c^. 
For this reason, we will examine the two-dimensional c^(0)-tan/3 plane for 
different values of Cg(0) and 0,^3(0). This naturally splits the four-dimensional 
space of initial values into two two-dimensional surfaces, greatly simplifying 
the analysis. 

20 



4.1 All m^ > 

Phenomenologically Allowed Regions and the Mass Spectrum: 

We first present our analysis subject to the following additional condition. 

• With the exception of ttlI,^, all squark and slepton soft squared masses 
are constrained to be positive over the entire scaling range. 

To illustrate the procedure, pick an arbitrary point 

Cg(0) = 0.75, 0^3(0) = 0.75 (69) 

in the Cq{0)-Ci^^{0) plane. For these initial values, we scan over the c^(0)-tan/3 
plane, first imposing the positive squark/slepton squared mass condition and 
then analyzing each point relative to the constraints discussed in the previous 
section. The results are shown in Figure iQ The positive squared mass 
condition is satisfied everywhere in the depicted region. 

Figure 1(a) shows the regions where electroweak symmetry is and is not 
radiatively broken, indicated in yellow and white respectively. The yellow re- 



gion is defined as the locus of points where both inequalities (64) and (65 ) are 
satisfied, whereas in any white region either one or both of these inequalities is 
violated. Before analyzing the individual areas, let us recall the consequences 



of each inequality. As discussed in Subsection 3.8, (64) guarantees that one 



linear combination of Higgs fields has a negative squared mass. In this case. 



satisfying inequality (65) implies a stable electroweak breaking vacuum. If, 



however, (65) is violated, the potential energy is not bounded from below 



and no stable vacuum state exits. On the other hand, violating inequality 



(64) indicates that the origin of Higgs space is either a local minimum or a 



local maximum of the potential energy, depending on whether or not (65) is 
satisfied. 

Let us now discuss the individual regions. Anywhere in the yellow region 



both (64) and (65) are satisfied, leading to a stable electroweak breaking 
vacuum. Note that there are two separated areas where electroweak breaking 
does not occur. Our analysis shows that at any point in the upper white 



^Here, and throughout the remainder of this paper, the irregularity of the plotted lines 
in the c^(0)-tan/3 and Cq{0)-c,y^{0) planes reflects the granularity of the grid structure used 
in the calculation. Each Cp(0)-tan/3 figure roughly represents 120^ data points, whereas 
there are roughly 200^ points in each Cq{0)-c,^^{0) plot. 
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(a) 



(b) 




20 25 iO 



(c) 



(d) 




(e) 

r IgUre li The c,j(0)-tan/3 plane corresponding to the point Cq{0) = 0.75, c^g (0) = 0.75. The yellow and 
white regions of (a) indicate where electroweak symmetry is and is not broken respectively. The individual 
regions satisfying the present experimental bounds for squarks and sleptons, gauginos and Higgs fields are 
shown in (b),(c) and (d), while their intersection is presented in (e). The dark brown area of (e) is 
the phenomenologically allowed region where electroweak symmetry is broken and all experimental mass 
bounds are satisfied. We present our predictions for the sparticle and Higgs masses at point (P). 
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region it is the first inequality (64) that is violated, while (65) continues 
to be satisfied. This indicates a stable vacuum, but with vanishing Higgs 
VEVs. The transition between the yellow and upper white regions is defined 



by saturating inequality (64), that is 



4 



(ic,r+4) 



\c>.\' + cl) 



(70) 



It follows from this and expression (60) that the boundary between these 
regions corresponds to the vanishing of M| in (61), that is. 



12 2 1 



^l-sin\2f3) 



Ch-c_ 



2|c„ 







(71) 



plotted as a function of tan (3 and c^(0). Below this boundary M^ is positive, 
indicating electroweak symmetry breaking vacua. At and above this line, 
however, M'^ vanishes, implying that electroweak symmetry is unbroken. 
Similarly, the lower right white region shown in Figure 1(a) also violates 



constraint (64) while satisfying (65). Hence, the above analysis applies here 



as well. For completeness, we point out that, beyond the boundaries shown 
in Figure 1(a), there is a transition of this lower right region to an area where 



both inequalities (64) and (65) are violated. In this regime, there are no stable 
vacua. 

Figures 1(b), (c) and (d) indicate where our calculated masses of the 
squarks, sleptons, Higgs and gauginos respectively exceed the experimen- 
tal lower bounds presented in Table 1. Finally, Figure 1(e) superimposes all 
of these with the area of electroweak symmetry breaking, the dark brown re- 
gion representing their intersection. Any point in this region has broken elec- 
troweak symmetry and a mass spectrum satisfying all experimental bounds. 
As an example, consider the point (P) indicated in this region. Our calcu- 
lated values for the squark, slepton, Higgs and gaugino masses are presented 
in Table 2. Note that, as stated, their values all exceed the experimental 
bounds. 



The above analysis was carried out for the arbitrarily chosen point (69) 
in the Cq{Q)-Cy^{^) plane. We emphasize that although this point has a 
non- vanishing region in the c^(0)-tan/3 plane satisfying all phenomenolog- 
ical bounds, this need not be the case for other points. To explore this, we 
now scan over the entire Cq{Q)-Cy^{^) plane. At each point, we analyze the as- 
sociated c^(0)-tan/3 plane and see if an allowed region exists. The results are 
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iable 2; The predicted spectrum at point (P) in Figure 1(e). The tilde denotes the superpartner of the 
respective particle. The superpartners of left-handed fields are depicted by an upper case label whereas 
the lower case is used for right-handed fields. The considerable mixing between the third family left- and 
right-handed scalar fields is incorporated into these results. 



shown in Figure [2|p| The white region indicates points whose corresponding 
c^(0)-tan/3 plane contains no locus of electroweak symmetry breaking. The 
yellow area represents points whose c^(0)-tan/3 plane has a region where 
electroweak symmetry is broken. Finally, each point in the blue area has 
a phenomenologically allowed region in its corresponding c^(0)-tan/3 plane 



satisfying the squark/slepton positive squared mass condition. Point (69) 



analyzed above is indicated by (A) in the diagram. It is of interest to see 
how the results change as we move to different phenomenologically allowed 
points in the Cq{0)-Ci,^{0) plane. For example, consider point (B) shown in 
Figure |2} This has the values 



Cg{0) = 1.4 , c,(0) = 1.2 



(72) 



For this point, the regions of the c^(0)-tan/3 plane corresponding to the 
different constraints, as well as their intersection, are shown in Figure |3] 



^ Note that the blue shaded allowed region can get infinitesmally close to, but not 
touch, both the horizontal and vertical axes. The somewhat irregular boundary lines 
reflect both the complexity of solving many RGEs, as well as the numerical limitations 
of our calculation. These comments apply to all figures of the Cq{Q)-Ci,^{0) plane in this 
paper. 
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C (0) 



r igUre L'. A plot of the Cq(0)-c^3(0) plane showing physically relevant areas. The yellow and white 
indicate points whose corresponding c^ (O)-tan fi plane does and docs not contain a region of electrowcak 
symmetry breaking respectively. Within the yellow area, the blue shading contains all points whose c^(0)- 
tan/3 plane has a non-vanishing region satisfying all experimental sparticle and Higgs bounds and for 
which all soft susy breaking masses remain positive over the entire scaling range. (A) and (B) indicate 
the two points analyzed in detail in the text. 

The positive squared mass condition is satisfied everywhere in the depicted 
regime. 



In the yellow region both (64) and (65) are satisfied, leading to stable elec 



troweak breaking vacua. There are two separated areas where electroweak 
breaking does not occur. As occurred for point (A), anywhere in the upper 



white region the first inequality (64) is violated, while (65) continues to be 



satisfied. This indicates stable vacua, but with vanishing Higgs VEVs. As 
discussed above, the boundary between the yellow and upper white regions 



corresponds to the vanishing of M^ in (61). Unlike the analysis of point 



(A), however, the lower right white region shown in Figure 3 violates both 



constraints (64) and (65). Hence, the origin of Higgs space is a local maxi- 



mum and the potential energy is unbounded from below. There are no stable 
vacua in this regime. 

The regions where the squarks/sleptons, gauginos and Higgs exceed their 
experimental lower bounds are depicted in the indicated colors. Any point in 
the intersection area, shown in dark brown, has broken electroweak symmetry 
and a mass spectrum satisfying all experimental bounds. As an example, 
consider the point (Q) indicated in this region. Our calculated values for the 
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r igUre o: The C;j(0)-tan/3 plane corresponding to the point Cq(0) = 1.4, c^3(0) = 1.2. The yellow 
and white regions indicate where electroweak symmetry is and is not broken respectively. The individual 
regions satisfying the present experimental bounds for squarks and sleptons, gauginos and Higgs fields 
are shown in the indicated colors. The dark brown area is their mutual intersection where electroweak 
symmetry is broken and all experimental mass bounds are satisfied. We present our predictions for the 
sparticle and Higgs masses at point (Q). 

squark, slepton, Higgs and gaugino masses are presented in Table 3. Note 
that, as stated, their values all exceed the experimental bounds. 

The i?-L/Electroweak Hierarchy: 

We have determined the subspace of the Cq(0)-Cy3(0) plane for which each 
point has a region in the corresponding c^(0)-tan/3 plane satisfying 1) the 
positive squark/slepton squared mass condition with 2) broken electroweak 
symmetry and 3) phenomenologically acceptable squark, slepton, Higgs and 
gaugino masses. Given such a point in the Cq(0)-Ci,3(0) plane and choosing a 
point in the acceptable region in the c^(0)-tan/3 plane, we now analyze the 
following question: What is the S-L/electroweak hierarchy for these initial 
values? 

An expression for the i?-L/electroweak hierarchy in terms of the q coef- 
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iable O! The predicted spectrum at point (Q) in Figure 3. The tilde denotes the superpartner of the 
respective particle. The superpartners of left-handed fields are depicted by an upper case label whereas 
the lower case is used for right-handed fields. The mixing between the third family left- and right-handed 
scalar fields is incorporated. 



ficients and tan/3 was given in (69). We repeat it here for convenience. 
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A/l-sin2(2/3) 



^H 



^H ^ I '^A' I 



1/2 



(73) 



For the specific point chosen in the initial Cg(0), 01,3(0), c^(0), tan /3 param- 
eter space, one can scale all quantities down to the electroweak scale and 
evaluate the hierarchy using (73). As a concrete example, consider point 



(A) in the Cq{0)-Ciy.^{0) plane of Figure 2. The corresponding regions of the 
c^(0)-tan/3 plane were superimposed in Figure 1(e) and are presented again 
in Figure 4(a). The allowed region is the dark brown area. For (A) given 



in (70), the 5-L/electroweak hierarchy is evaluated for each point in this 
allowed region and plotted in Figure 4(b). We find that the hierarchy takes 
values of 6.30-6.36 along the lower boundary of the allowed region. Note 
that below this boundary at least one of the gaugino or Higgs masses vi- 
olates their experimental bound. Hence, the lower values of the hierarchy 
are determined from the experimental data. On the other hand, as one ap- 
proaches the boundary with the upper white region, the hierarchy becomes 



infinitely large. To understand this, recall from (71) that this boundary is 
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r igUre 4; Plot (a) shows the c^(0)-tan/3 plane corresponding to point (A) in Figure 2 with the 
phenonienologically allowed region indicated in dark brown. The mass spectrum at (P) was presented in 
Table 2. A plot of the hierarchy Mg^i^/Mz over the allowed region is given in (b). Graph (c) shows the 
hierarchy as a function of c^(0) along the tan/3 = 18 line passing through (P). 
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determined by the vanishing of M| in ( 62 ) , that is 



'"^ '^' 4-4-2|c,r = 0. (74) 



Vl-sin2(2/3) 



Hence, at any point on this boundary the denominator in (73) vanishes and 



M, 



oo . (75) 



It follows that within the phenomenologically acceptable region, any value 
of the B-L hierarchy in the range 6.30 < Mab_lIMz < C)0 can he attained. 
Another way to analyze this data is to pick a specific point in the al- 



lowed region and to compute (73) as a function of c^(0) along the fixed tan/3 
line passing through it. For concreteness, choose the point (P) for which 
we calculated the mass spectrum in Table 2. This is shown in Figure 4(a) 
along with the dotted line tan/3 = 18 intersecting it. The 5-L/electroweak 
hierarchy along this line is plotted in Figure 4(c). Note that this begins at 
Mab_i^/Mz = 6.35 at the experimentally determined lower boundary, rises 
slowly to Mab_,^/Mz ~ 20 across most of the region, and then rapidly di- 
verges to infinity as one approaches the upper boundary. Approaching both 
the lower and, especially, the upper boundary requires fine-tuning of c^(0). 
For "typical" values of c^(0), the hierarchy is naturally in the range 

10 < ^^%^ < 20 . (76) 

^ Mz 

As a second example, consider point (B) in the Cq{Q)-Cy.^{Q) plane of Figure 
2. The corresponding regions of the c^(0)-tan/3 plane were superimposed in 
Figure 3 and presented again in Figure 5(a). The allowed region is the 



dark brown area. For (B) given in (73), the 5-L/electroweak hierarchy is 



evaluated for each point in this allowed region and plotted in Figure 5(b). We 
find that the hierarchy takes values of 10.00-10.21 along the lower boundary 
of the allowed region, below which at least one of the gaugino or Higgs masses 
violates their experimental bound. Again, as one approaches the boundary 
with the upper white region, the hierarchy becomes infinitely large. It follows 
that within the phenomenologically acceptable region any value of the B-L 
hierarchy in the range 10 < Mab_,^/Mz < oo can he attained. 

Another way to analyze this data is to pick a specific point in the al- 



lowed region and to compute (73) as a function of c^(0) along the fixed tan/3 
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r lgUr6 O! Plot (a) shows the c^(0)-tan/3 piano corresponding to point (B) in Figure 2 with the 
phenoinenologically allowed region indicated in dark brown. The mass spectrum at (Q) was presented in 
Table 3. A plot of the hierarchy Mg^i^/Mz over the allowed region is given in (b). Graph (c) shows the 
hierarchy as a function of c^(0) along the tan/3 = 12 line passing through (Q), 
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line passing through it. For concreteness, choose the point (Q) for which 
we calculated the mass spectrum in Table 3. This is shown in Figure 5(a) 
along with the dotted line tan/3 = 12 intersecting it. The 5-L/electroweak 
hierarchy along this line is plotted in Figure 5(c). Note that this begins at 
^Ab-lI^z = 10.15 at the experimentally determined lower boundary, rises 
slowly to Mab_,^/Mz ~ 30 across most of the region, and then rapidly di- 
verges to infinity as one approaches the upper boundary. For "typical" values 
of c^(0) not fine-tuned near either boundary, the hierarchy is naturally in the 
range 

15 < ^^%^ < 30 . (77) 

^ Mz 

4.2 m^^ < 

The Potential Energy for m^,^ < and mg^ < 0: 



For the choice of parameters in (55), all sleptons have positive soft squared 



masses with the exception of the third family right-handed sneutrino, for 

which m1^ < 0. As noted in Subsection 3.8, imposing positivity on the 

effective masses of the left-handed squarks at the B-L breaking VEV (z/3), 

that is, 

1 

does not require that rn^, be positive. In general, one or more of these 



(mjj=m^, + -^4^(^/3)' >0, (78) 



soft squared masses can be negative. Despite our assumption in (20), (35) 
that the initial squark masses are universal, the effect of the large third 
family up- Yukawa coupling in the RGEs is to break this degeneracy, driving 
ttlq^ negative more quickly than the first and second family squark masses. 
Therefore, for simplicity, we explore the possibility that only the third family 
left-handed squark soft mass becomes negative, toq^ < 0, as it is scaled down 
to electroweak energy-momenta. 

The electroweak phase transition breaks the left-handed SU{2)l doublet 
Qs into its up- and down- quark components f/3 and D^ respectively. The 
leading order contribution of the Higgs VEVs to their mass splits the degen- 
eracy between these two fields, destabilizing the potential most strongly in 
the D^ direction. For this reason, the relevant Lagrangian for analyzing this 
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vacuum can be restricted to 

4 



>C = iV^^^us 



FB-LiJ.uFg_j^ + \'Dj:,^^D-i\ 



Y fiu^Y 



Fr 



-^Fsu{2)i,uF^lri2) - 4^5[/(3Vi.^|^(3) - '^^(^^3,^3) 



where 



V 



Will — d^ — igB-L^B-Lfi 



(79) 



(80) 



^D3M 



r, .Qb-l , -Sy . ■ . ■ . 

Oil - i^-^B-Lf, - ^-^^Yfi - ig2AsU(2)t, - «5'3^5J/(3)a 



and 



V{U3,D3 



2 I 1 2 I 2 I T^ 



1 



"1/31 



31^ +%^(|-3r+ 31/^3 



|2\2 



(81) 



1/fi'y 5'2 
+ -(— + — + 
2^36 4 3 



5'3mt-, |4 



The first two terms in the potential are the soft supersymmetry breaking 
masses in (12), while the remaining terms are supersymmetric and arise 
from Db-l, Ffy in ([lO]), (jO)) and Dsu(2)l, F>su{3)c respectively. Using Xa^ ~ 
5 X 10^^, a hierarchy with (IF^) ^ {1^3) and assuming jm^gl is of order |r«,^3|, 
terms proportional to the Higgs VEVs are small and are ignored in (81). For 
simplicity, we henceforth drop the small (y'^_^/9 + (7y/36 piece of the D-term 
contribution. 

If both m^g < 0, mjj^ < at the electroweak scale, then the potential is 
unstable at the origin of field space and has two other local extrema at 



(^3)^ 



m 



vi 



2 ' 
Tb-l 



(D-s) = 



and 



(z/3) = 0, {D^y 



m 



Ds 



9l/^ + ^1/3 



respectively. Using these, potential (|81|) can be rewritten as 

^(^^3,^3 



9I-L 



(k3| 



/ \2\2 I 9b-L I |2| n |2 



I ^1/4 + ^i/3 



(lA 



3 



(82) 



(83) 



(84) 
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Let us analyze these two extrema. Both have positive masses in their ra- 
dial directions. At the sneutrino vacuum (82), the mass squared in the D^ 
direction is given by 

|2 



^DzWi) 



9I-L 



>^y 



;f + f)(^3)^ 



\m,. 



IniDs 



(85) 



whereas at the D3 vacuum (83), the mass squared in the z/3 direction is 



771 



^z 



(Ds) 



gl^L 



{D.f-gVd^z? 



F^D3 



gl^L 



'^9l/^ + 9l 



— \m 



Ki\ 



(86) 



Note that either (85) or (86) can be positive, but not both. To be consistent 



with the hierarchy solution, we want (82) to be a stable minimum. Hence, 



we demand fn'j~,J(^i,^) > or, equivalently, that 



m,. 



> 3|m/3, 



(87) 



We will impose (|87|) as an additional condition for the remainder of this 

< and, hence, the 



subsection. It then follows from (86) that m 



1^3 



\{D3) 



D^ extremum (83) is a saddle point. As a consistency check, note that 



y|(i^3) < V\(^D3) if and only if 



^Ll(^3)' > ( 
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9l 



){D,Y 



or, equivalently, 



m,. 



> k«D3 



9I-L 



\l/2 



(88) 



(89) 



'^9l/^ + 9V 
This follows immediately from constraint ( [87) ). 

Finally, note that the potential descends monotonically along a path C 



from the saddle point at (|83|) to the absolute minimum at (|84|). Solving the 

dV 
dDa 



equation, this curve is found to be 
((^3)^ 



D 



3\C 



\l^3 



9%-L 



= )) 



1/2 



(90) 



Note that it begins at {D3) for u^ 
intersects the D3 = axis at |z/3o| 



and continues until it tangentially 
I I xZ/s). From here, the path 



^\rnn,\ 



continues down this axis to the stable minimum at (82). We conclude that 



at the electroweak scale the absolute minimum of potential (81) occurs at 



the sneutrino vacuum given in (82). 
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Phenomenologically Allowed Regions and the Mass Spectrum: 

In this subsection, we analyze our results subject to the following additional 
conditions. 

• The third family left-handed down-squark soft mass squared will be con- 
strained to be negative, that is, mfj^ < 0. All other squark and slepton 
soft squared masses are positive over the entire scaling range, with the 
exception of m^^ . 

• To ensure that the B-L breaking VEV is the absolute minimum, we 
impose condition (87), 

|m,^3p > 3|m£)3p , (91) 

at the electroweak scale. 

We will refer to these two conditions collectively as the ?7i|,^ < mass con- 
dition. 

As discussed in the previous subsection, we proceed by scanning over the 
entire Cq{0)-c,^^{0) plane, at each point analyzing the associated c^(0)-tan^ 
plane to see if an allowed region exists. The results are shown in Figure 6. 
As in Figure 2, the white region indicates points whose corresponding c^(0)- 
tan/3 plane contains no locus of electroweak symmetry breaking, whereas 
the yellow area represents points whose c^(0)-tan/3 plane has a region where 
electroweak symmetry is broken. Finally, each point in the red area has 
a phenomenologically allowed region in its corresponding c^(0)-tan/3 plane 
satisfying the mj^^ < mass condition. Note that this is distinct from the 
blue region in Figure 2, where all squark/slepton mass squares are positive. 
Let us analyze the properties of an arbitrary point in the red area. For 
example, consider point (C) shown in Figure 6. This has the values 

Cg(0) = 1.0 , c,(0) = 1.1 . (92) 

For this point, the regions of the c^(0)-tan/3 plane corresponding to the 
different constraints, as well as their intersection, are shown in Figure 7. The 
m'jj^ < mass condition is satisfied everywhere in the depicted regime. 



In the yellow region both (64) and (65) are satisfied, leading to stable 
electroweak breaking vacua. There are two separated areas where electroweak 
breaking does not occur. As for point (A) in Figure 2, anywhere in the upper 
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0.25 03 0.75 1,„, 1.25 15 1.75 

r lgU.r6 u! A plot of the Cq(0)-Ci/.^(0) plane showing physically relevant areas. The yellow and white 
indicate points whose corresponding c^ (O)-tan plane does and does not contain a region of electroweak 
symmetry breaking respectively. Within the yellow area, the red shading contains all points whose c^(0)- 
tan/9 plane has a non-vanishing region satisfying all experimental sparticlc and Higgs bounds and for 



and lower right white regions the first inequahty (64) is violated, while (65) 
continues to be satisfied. This indicates stable vacua, but with vanishing 
Higgs VEVs. It follows that the boundary between the yellow and white 



regions corresponds to M| in (61) becoming zero. The regions where the 
squarks/sleptons, gauginos and Higgs exceed their experimental lower bounds 
are depicted in the indicated colors. Any point in the intersection area, shown 
in dark brown, has broken electroweak symmetry and an acceptable mass 
spectrum. As an example, consider the point (R) indicated in this region. 
Our calculated values for the squark, slepton, Higgs and gaugino masses are 
presented in Table 4. 



The 5-L/Electroweak Hierarchy: 

We have determined the subspace of the Cg{0)-c,^^{0) plane for which each 
point has a region in the corresponding c^(0)-tan/3 plane satisfying 1) the 



m 



2)jj < mass condition with 2) broken electroweak symmetry and 3) phe- 
nomenologically acceptable squark, slepton, Higgs and gaugino masses. Given 
such a point in the Cq{0)-c^.^{0) plane and choosing a point in the acceptable 
region in the c^(0)-tan/3 plane, we now analyze the i?-L/electroweak hierar- 
chy for these initial values. 
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r igUre i'. The c^(0)-tan/3 plane corresponding to the point Cq(0) = 1.0,c^3(0) = 1.1. The yellow 
and white regions indicate where electroweak symmetry is and is not broken respectively. The individual 
regions satisfying the present experimental bounds for squarks and sleptons, gauginos and Higgs fields 
are shown in the indicated colors. The dark brown area is their mutual intersection where electroweak 
symmetry is broken and all experimental mass bounds are satisfied. We present our predictions for the 
sparticle and Higgs masses at point (R). 
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iable 4; The predicted spectrum at point (R) in Figure 7. The tilde denotes the superpartner of the 
respective particle. The superpartners of left-handed fields are depicted by an upper case label whereas 
the lower case is used for right-handed fields. The mixing between the third family left- and right-handed 
scalar fields is incorporated. 
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An expression for this hierarchy in terms of the Cj coefficients and tan /3 



was given in (73 ). For the specific point chosen in the initial Cq{0) , Ci,^{Q) , c^(0), 
tan/3 parameter space, one can scale all quantities down to the electroweak 
scale and use this expression to evaluate the hierarchy. As a concrete ex- 
ample, consider point (C) in the Cq{0)-c^.^{0) plane of Figure 6. The corre- 
sponding regions of the c^(0)-tan/3 plane were superimposed in Figure 7 and 
are presented again in Figure 8(a). The allowed region is the dark brown 



area. For (C) given in (92), the 5-L/electroweak hierarchy is evaluated for 
each point in this allowed region and plotted in Figure 8(b). We find that 
the hierarchy takes values of 8.99-9.06 along the lower boundary of the al- 
lowed region. Note that below this boundary at least one of the gaugino 
or Higgs masses violates their experimental bound. Hence, the lower values 
of the hierarchy are determined from the experimental data. On the other 
hand, as one approaches the boundary with the upper white region, the hi- 
erarchy becomes infinitely large. As discussed in the previous subsection. 



this is explained by the vanishing of M^ in (61). Hence, at any point on 



this boundary the denominator in (73) vanishes and MAg_^/Mz — > oo. It 



follows that within the phenomenologically acceptable region, any value of 
the B-L hierarchy in the range 8.99 < Mab_^IMz < C)0 can he attained. 
Another way to analyze this data is to pick a specific point in the al- 



lowed region and to compute (73) as a function of c^(0) along the fixed tan/3 
line passing through it. For concreteness, choose the point (R) for which 
we calculated the mass spectrum in Table 4. This is shown in Figure 8(a) 
along with the dotted line tan/3 = 14 intersecting it. The i?-L/electroweak 
hierarchy along this line is plotted in Figure 8(c). Note that this begins at 
^Ab-lI-^z = 9.0 at the experimentally determined lower boundary, rises 
slowly to Mab_,^/Mz ~ 40 across most of the region, and then rapidly di- 
verges to infinity as one approaches the upper boundary. Approaching both 
the lower and, especially, the upper boundary requires fine-tuning of c^(0). 
For "typical" values of c^(0), the hierarchy is naturally in the range 

15 < ^^%^ < 40 . (93) 

^ My. 



"Mixed" rri^ > and mj^^ < Mass Conditions: 

It is of interest to superimpose the blue region in Figure 2, satisfying the m^ > 
mass condition, with the red region of Figure 6, defined by the m'jy, < 
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r IgUre o: Plot (a) shows the Cf,(0)-tan/3 plane corresponding to point (C) in Figure 6 with the 
phenomenologically allowed region indicated in dark brown. The mass spectrum at (R) was presented in 
Table 4. A plot of the hierarchy Mg_ij/Mz over the allowed region is given in (b). Graph (c) shows the 
hierarchy as a function of c^(0) along the tan/3 = 14 line passing through (R). 
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r igUre 9; A plot of the Cq (0)-c^3 (O) plane showing both the blue and red regions presented in Figures 2 
and 6 respectively. They have a non-vanishing intersection, indicated in purple. Any point in this overlap 
has an allowed region in the Cp(0)-tan/3 plane that is divided into two areas-one with rr? > and the 
second with m,|j < 0. (D) indicates a point in this overlap region analyzed in detail in the text. 



constraint. This is shown in Figure 9. Note that there is a non-vanishing 
intersection between these two areas. This is comprised of points in the Cg(0)- 
01,3(0) plane whose phenomenologically allowed regions in the corresponding 
c^(0)-tan/3 plane are each divided into two regimes-one satisfying the m^ > 
mass condition and the other the vn\)^ < constraint. As a specific example, 
consider the point (D) shown in Figure 9. This has the values 



c,(0) = 1.0 , c,(0) = 0.9 



(94) 



For this point, the areas of the c^(0)-tan/3 plane corresponding to the dif- 
ferent constraints, as well as their intersection, are shown in Figure 10. The 
regions where the squarks/sleptons, gauginos and Higgs exceed their exper- 
imental lower bounds are depicted in the indicated colors. Any point in the 
intersection area, shown in dark brown, has broken electroweak symmetry 
and an acceptable mass spectrum. Importantly, however, note the dotted 
line dividing this plane. We find that the m"^ > mass condition is satisfied 
everywhere to the left of this line, whereas the m'jj, < constraint holds at 
all points to the right-consistent with (D) being a point in the intersection 
of the blue and red regions. The dotted line is vertical since, to leading order 
the D3 mass squared, although a function of tan/3, is independent of c^(0). 
The sparticle and Higgs mass spectrum for point (S) in the allowed region is 
presented in Table 5. 
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Figure 10: The Cp(0)-tan/3 plane corresponding to the point Cq{0) = 1. 0,0^3(0) = 0.9. The yellow 
and white regions indicate where electroweak symmetry is and is not broken respectively. The individual 
regions satisfying the present experimental bounds for squarks and sleptons, gauginos and Higgs fields 
are shown in the indicated colors. The dark brown area is their mutual intersection where electroweak 
symmetry is broken and all experimental mass bounds are satisfied. We present our predictions for the 
sparticle and Higgs masses at point (S). The dotted line passing to the right of (S) separates the m? > 
region, to the left of this line, from the area where m|j < 0, to the right. 
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iable 5; The predicted spectrum at point (S) in Figure 10. The tilde denotes the superpartner of the 
respective particle. The superpartners of left-handed fields are depicted by an upper case label whereas 
the lower case is used for right-handed fields. The mixing between the third family left- and right-handed 
scalar fields is incorporated. 



40 



4.3 mi <0 

The Initial Conditions and Potential Energy for m^^ < and ml^ < 0: 

As discussed in Subsection 3.8, to guarantee that the B-L vacuum is a stable 
local minimum, we impose the constraint that the effective squared masses 
of all squarks and sleptons evaluated at (z/3) are positive over the entire 



scaling range. Similarly to the left-handed squark mass condition (77), im- 
posing positivity on the effective right-handed down slepton masses at the 
B-L breaking VEV (z/3), that is, 

(ml) = ml + ^gl{us)' > , (95) 

does not require that m-g, be positive. In general, one or more of these soft 



squared masses can be negative. Recall from (21 ),(36) that we have assumed 
that all left-handed and right-handed down sleptons have a universal initial 
mass. This is similar to the initial condition on squark masses. Unlike the 
squarks, however, the down- Yukawa couplings of sleptons are all too small 
to greatly effect the RGE running of their soft masses. It follows that, at 
a low scale, the three families of right-handed down sleptons mass squares 
tend to be all positive or all negative. Splitting this degeneracy, for example, 
to drive only ml,^ < 0, requires considerable fine-tuning. Therefore, if one 
wishes to consider the case where only the third family squared mass turns 
negative, it is necessary to alter the initial slepton mass conditions given in 
Section 3. This is easily accomplished as follows. 

As discussed in Subsection 3.3, the boundary conditions for the RGEs of 
the Higgs, squarks and sleptons squared masses are greatly simplified if one 
chooses the initial soft masses so that both 5(0) = and iSq(O) = 0, with S 



and Sq given in (19) and (23) respectively. Hence, in this paper we always 
choose the initial parameters to satisfy these two conditions. However, the 
specific choices made in Subsection 3.3 were overly constraining, since they 
imposed unification of all three families of squarks and sleptons, whereas the 
unification of each family separately is sufficient. In particular, condition 



(|21j) sets 

mLX0f = m,^{0f (96) 

for all i,j = 1,2,3. This leads to the difficulty discussed above. However, 
this constraint can easily be weakened. The simplest example is to take 

mi,,,(0)2 = m,,M' , ^lM' = ^eM' (97) 
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which clearly continues to solve both 5(0) = and 5q(0) = 0. Expression 



(36) then generalizes to 

rriL^^iO) = me^2(0) = Ce,^^{0)M , tul^ = me,{0) = Ces{0)M . (98) 



In terms of these parameters, (45) becomes 

5^(0) = {1 + 2C^ 
where 



C 



^;^i,2 



(0) 



A 






(99) 



(100) 



CuM ' ^^ c^3(0) ' ^^^ 0^.3(0) ■ 
To stay as close as possible to our previous analysis, we continue to use the 



values 



A = VQ, C = 9.12 



introduced in (51) and (53) respectively. In addition, let us choose 

A3 = Vs , 



(101) 



(102) 



thus minimally changing the value of (47) from 5 to 5.1. It follows that 



equations (48), (54), and the conclusions thereof for U{1)b-l breaking, do 



not change substantially. Similarly, equation (52) for i = 1, 2 is minimally 
altered to 

m,^,{tB-L? = {{V6f - 5.1) m,3(0)2 = 0.9 m,^{Of . (103) 

However, we now find that 

me,{tB-Lf = {{Vsf - 5.1) m,,{Of = -2.1 m,^{Of . (104) 

That is, splitting the slepton coefficient into A = y/E and A^ = a/3 allows the 
mass squares of the first two families to remain positive while constraining 



TBg < 0, as desired. Henceforth, (55) is replaced by 



Cei,2(0) = Vg 0,3(0) , Ce3(0) = 73 0,3(0) , c,,M = 9.12 0,3(0) . (105) 
Despite these changes in the initial conditions, Cg(0), 0,3(0), c^(0) and tan/3 



in (63) remain the four independent parameters of our analysis. 



The new set of initial parameters just discussed allows for the possibility 
that, at the electroweak scale, all soft squared masses are positive with the 
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exception of ml < and m^, < 0. The relevant potential for discussing the 



vacuum of z/3 and 63 is given by 

T// N 2 I |2 I 2 I |2 I 9b-L i\ I 

v/(«^3,e3) =m^Ji/3| +mgje3| H ;^— (ksl 



leslY + fhr- (106) 



631-1 ■ 2 

The first two terms in the potential are the soft supersymmetry breaking 



mass terms in (12), while the third and fourth terms are supersymmetric 



and arise from the Db-l and Dy in (10) and ^ respectively. Contributions 
to ( 106 ) from the relevant Yukawa couplings in ^ are suppressed, since X^^ 



and Ae3 are of order 10^^° and 10^^ respectively. Hence, we ignore them. 

If both ml,^ < 0, TTigg < at the electroweak scale, then the potential is 
unstable at the origin of field space and has two other local extrema at 

2 

Til 

(^3)^ = -V^, (e3) = 



2 ' 
9b-l 



(107) 



and 



(^^3) = 0, (63)' 



m 



63 



9%-L + 9y 



(108) 



respectively. Using these, potential (|106|) can be rewritten as 
^(^^3,63) 



9I-L 



(ksl 



/ \2\2 I 2 I |2| |2 

(^^3) ) +9b-lW3\ m 



9l^L + 9y 



2\2 



(Ie3|^-(e3)^) 



(109) 



Let us analyze these two extrema. Both have positive masses in their 



radial directions. At the sneutrino vacuum (107), the mass squared in the 63 
direction is given by 



m^ 



K'^s) 



9l-L(yzf - {9l^L + 9l){(i-i) 



;iio) 



whereas at the stau vacuum (|108|), the mass squared in the v^ direction is 

(111) 



'i, 



m.,„|(e3) 



9B-L{e3) -9B-LKV?>y 



m, 



63 1 



1 + 



9I 
9I-L 



m 



2^3 1 



Note that either ( 110 ) or ( 111 ) can be positive, but not both. To be consistent 



with the hierarchy solution, we want (|107|) to be a stable minimum. Hence, 

i^ I 

'63 1 

;ii2) 



we demand m1^ \ (^,3) > or, equivalently, that 



\my^\^ > Ime^l"^ 
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We will impose (112) as an additional condition for the remainder of this 

< and, hence, the 



subsection. It then follows from (111) that m 



^3 



I(e3> 



stau extremum (108) is a saddle point. As a consistency check, note that 



1^1 (j^g) < ^1(63) if and only if 



9l-L{^^)'>{9l^L + 9Y){e,Y 



or, equivalently, 



771 



V3\ 



> \m 



63 1 



:i + 



9y - 
9%-l' 



-1/2 



;ii3) 



;ii4) 



This follows immediately from constraint (112). Finally, note that the po- 



to the absolute minimum at (|107|). Solving the |^ 



tential descends monotonically along a path C from the saddle point at (|108|) 
to the absolut 
is found to be 



leslc = ((es)^ - I 

Note that it begins at (e^) for u^ 
intersects the 63 = axis at |z/3o| = 



^3r{i + 



9y - 
9I-L 



-Ul/2 



equation, this curve 
(115) 



and continues until it tangentially 
I's). From here, the path continues 



Im^gl 



down this axis to the stable minimum at (95). We conclude that at the 



electroweak scale the absolute minimum of potential (106) occurs at the 



sneutrino vacuum given in (107). 



Phenomenologically Allowed Regions and the Mass Spectrum: 

In this subsection, we analyze our results subject to the following additional 
conditions. 

• The third family right-handed slepton soft mass squared will be con- 
strained to be negative, that is, m^^ < 0. All other squark and slepton 
soft squared masses are positive over the entire scaling range, with the 
exception of m^^ . 

• To ensure that the B-L breaking VEV is the absolute minimum, we 



impose condition (112) 



Imj^gl^ > Imegl^ 



(116) 



at the electroweak scale. 
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We will refer to these two conditions collectively as the ml,^ < mass condi- 
tion. 

As discussed in previous subsections, we proceed by scanning over the 
entire Cq{0)-c^^{0) plane, at each point analyzing the associated c^(0)-tan/3 
plane to see if an allowed region exists. The results are shown in Figure 11. 
As in Figures 2 and 6, the white region indicates points whose correspond- 
ing c^(0)-tan/3 plane contains no locus of electroweak symmetry breaking, 
whereas the yellow area represents points whose c^(0)-tan/3 plane has a re- 
gion where electroweak symmetry is broken. Finally, each point in the green 
area has a phenomenologically allowed region in its corresponding c^(0)-tan/3 
plane satisfying the ml,^ < mass condition. Since some of the initial param- 
eters are now different to allow for a negative stau squared mass, this green 
region cannot be superimposed with the blue and red regions discussed pre- 
viously. Let us analyze the properties of an arbitrary point in the green area. 
For example, consider point (E) shown in Figure 11. This has the values 

c,(0) = 1.1 , c,(0) = 0.5 . (117) 

For this point, the regions of the c^(0)-tan/3 plane corresponding to the 
different constraints, as well as their intersection, are shown in Figure 12. 
The r/igg < mass condition is satisfied everywhere in the depicted regime. 



In the yellow region both inequalities (64) and (65) are satisfied, leading 
to stable electroweak breaking vacua. There are two separated areas where 
electroweak breaking does not occur. As before, anywhere in the upper white 



region the first inequality (64) is violated, while (65) continues to be satisfied. 
This indicates stable vacua, but with vanishing Higgs VEVs. It follows that 
the boundary between the yellow and upper white regions corresponds to the 



vanishing of M| in (61). However, as at point (B), for example, the lower 



right white region shown in Figure 12 violates both constraints (64) and (65). 
Hence, the origin of Higgs space is a local maximum and the potential energy 
is unbounded from below. There are no stable vacua in this regime. 

The regions where the squarks/sleptons, gauginos and Higgs exceed their 
experimental lower bounds are depicted in the indicated colors. Any point in 
the intersection area, shown in dark brown, has broken electroweak symmetry 
and an acceptable mass spectrum. As an example, consider the point (T) 
indicated in this region. Our calculated values for the squark, slepton, Higgs 
and gaugino masses are presented in Table 6. 
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0.25 0,5 0,75 






T.75 



r igUre 11; A plot of the Cq(0)-Ci/3(0) piano showing physically relevant areas. The yellow and white 
indicate points whose corresponding C(j (O)-tan fi plane does and does not contain a region of electroweak 
symmetry breaking respectively. Within the yellow area, the green shading contains all points whose 
Cji (0)-tan/3 plane has a non- vanishing region satisfying all experimental sparticle and Higgs bounds and 
for which m^ < 0. (E) indicate the point analyzed in detail in the text. 



1.7 
1.6 
1.5 
Cf, (0) 1.4 
13 
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r igUre \L'. The Cfi(0)-tan/3 plane corresponding to the point Cq(0) = 1.1,0^3(0) = 0.5. The yellow 
and white regions indicate where electroweak symmetry is and is not broken respectively. The individual 
regions satisfying the present experimental bounds for squarks and sleptons, gauginos and Higgs fields 
are shown in the indicated colors. The dark brown area is their mutual intersection where electroweak 
symmetry is broken and all experimental mass bounds are satisfied. We present our predictions for the 
sparticle and Higgs masses at point (T). 
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Particle 


Symbol 


Mass [GeV 


Particle 


Symbol 


Mass [GeV] 


S quarks 


^Ql,2 

w\br 

f(l) f(2) 
'-3 ' ''3 


1176 

1136, 1184 

932, 1057 

770, 986 


Higgs 




132 
640 
640 

645 


Sleptons 


^1,2 

n,2 
^3 '^3 


806 

768 

519, 606 


Neutralinos 




146 
290 
726 
734 


Charginos 


Al ) A2 


290, 734 


Gluinos 


9 


1066 


Z' Boson 


Ab-l,Ab-l 


1776, 1860 



iable 6; The spectrum at point (T) in Figure 12. The tilde denotes the superpartncr of the respective 
particle. The superpartners of left-handed fields are depicted by an upper case label whereas the lower 
case is used for the right-handed fields. The mixing between the third family left- and right-handed scalar 
fields is incorporated. 



The i?-L/Electroweak Hierarchy: 

We have determined the subspace of the Cq{0)-c,y.^{0) plane for which each 
point has a region in the corresponding c^(0)-tan/3 plane satisfying 1) the 



m 



63 



< mass condition with 2) broken electroweak symmetry and 3) phe- 
nomenologically acceptable squark, slepton, Higgs and gaugino masses. Given 
such a point in the Cq{0)-Ci„^{0) plane and choosing a point in the acceptable 
region in the c^(0)-tan/3 plane, one can analyze the 5-L/electroweak hier- 
archy for these initial values. The analysis proceeds exactly as in previous 
subsections, so we simply present the results. 

For point (E) in Figure 11, we have computed the hierarchy everywhere 
in the dark brown area of Figure 12. We find that this takes values of 
7.60-7.74 along the lower boundary of the allowed region. Note that be- 
low this boundary at least one of the gaugino or Higgs masses violates their 
experimental bound. Hence, the lower values of the hierarchy are deter- 
mined from the experimental data. On the other hand, as one approaches 
the boundary with the upper white region, the hierarchy becomes infinitely 
large for the reasons previously discussed. It follows that within the phe- 
nomeno logically acceptable region, any value of the B-L hierarchy in the 
range 7.60 < Ma^^lI^z < oo can he attained. 
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Another way to analyze this data is to pick a specific point in the allowed 
region and to compute (73) as a function of c^(0) along the fixed tan/3 line 



passing through it. For concreteness, choose the point (T) with tan/3 = 22 
for which we calculated the mass spectrum in Table 6. We find that the 
hierarchy begins at Mab_j^/Mz = 7.65 at the experimentally determined 
lower boundary, rises slowly to Mab_^/Mz ~ 30 across most of the region, 
and then rapidly diverges to infinity as one approaches the upper boundary. 
Approaching both the lower and, especially, the upper boundary requires 
fine-tuning of c^(0). For "typical" values of c^(0), the hierarchy is naturally 
in the range 

10<^^%^<30. (118) 

4.4 Summary 

We first note that the above classification of vacua using the sign of toq, 
and m\_ is complete. The only other squared masses are for right-handed 
squarks and left-handed sleptons, which enter the effective masses at the 
B-L breaking VEV {v^) as 

«) = < - 1^4(^3)' , K) = ml - \gl{u,f (119) 

and 

{ml:) = ml-\gl{v,)' (120) 

respectively. Since all of these effective masses must be positive to ensure 
that the vacuum is a stable minimum, it follows from the minus signs in each 
expression that m^ ,m^ , and m\, must all be positive. Therefore, all im? > 0, 
rriQ^ < 0, and ml. < in subsections 4.1, 4.2 and 4.3 respectively are the 
only possibilities. 

From the above analysis, several broad conclusions can be made. For the 
reasons discussed above, we limited our search to the four-dimensional space 



of parameters listed in ( 63 ) . By combining the results in the m'^ > , ttiq. < , 
and ml. < regimes, we can find the generic region of this parameter space 
for which one obtains a phenomenologically acceptable vacuum. The full 
range of allowed values for the Cg(0) and €^^{0) parameters were presented in 



Figures |9] and 11 From these, we observe a maximum range of 

< 0,(0) < 1.8, < 0^3(0) < 1.5 . (121) 
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Similarly, by examining the c^(0)-tan/3 plane over the allowed values of Cg(0) 
and 01,3(0), the range of phenomenologically allowed values is found to be 

0.8 < c^(0) < 1.75 , 8<tan;g<33. (122) 

To obtain this result, we computed the allowed regions for numerous points in 
the Cq{0)-Cu3{0) plane including, but not limited to, (A)-(E) presented in the 
text. Thus, even with our restrictive premises in Section 3, a phenomenolog- 
ically viable B-L MSSM vacuum exhibiting an acceptable hierarchy occurs 
for a reasonably wide space of initial parameters. Lifting some of the above 
constraints, such as allowing allmQ_ < instead of just ttIq < 0, will clearly 
allow a significant enlargement of the acceptable initial parameter space. 

5 Some (z/3) ^ Phenomenology 

The results presented in this paper allow one to compute any quantity in our 
B-L MSSM theory at any energy scale. In particular, we have shown that 
for a wide range of initial conditions there is a stable vacuum which breaks 
both B-L and electroweak symmetry with an acceptable sparticle and Higgs 
mass spectrum and 5-L/electroweak hierarchy. These are important nec- 
essary conditions on the theory, but are not sufficient to guarantee that it 
is phenomenologically viable. In this section, we explore two more impor- 
tant constraints arising from lepton number and baryon number violation 
respectively. 

5.1 Lepton Number Violation 

The most general superpotential invariant under gauge group SU(3)c x 
SU{2)l X U{1)y X U{1)b-l is presented in ^. Assuming a flavor diago- 
nal basis, the superpotential becomes 

3 
W = i^HH + Y^ {Xu,iQiHui + Xd,iQiHd + K^^LiHui + Xe,iL,Hei) . (123) 

Recall that since U{1)b-l contains matter parity, the dangerous lepton and 
baryon number violating terms in ([T]) are forbidden. Note, however, that 
these results are only valid at high scales where the gauge symmetry, in 
particular U{1)b-l, is exact. At low energy-momentum the gauged B-L 
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symmetry is spontaneously broken, potentially allowing these operators to 
"grow back". This can be analyzed by expanding the third family right- 
handed sneutrino around its VEV, that is, let z/3 = (z/3) + z/3. Note that 

fiHH + X,,LsHu3 = fiH{H + e3Ls) + ... , (124) 

where 

e, = xj-^. (125) 

/i 

This motivates performing a rotation of the down Higgs and third family 
lepton doublet superfields given, to leading order, by 

H' = H + esLs, L'^ = L3-e^H. (126) 

Written in terms of these new superfields, and then dropping the ' for sim- 
plicity, the superpotential becomes 

3 3 

W = W + es^Ae.iLsLiCi + e3J2^d,iL3Q^di , (127) 

i=l i=l 



where W is given in (123). As expected, the lepton number violating terms 
of the form 

LsL^e, , LsQidi (128) 

have grown back. Note, however, that the baryon violating terms Uidjdk have 
not been regenerated by the right-handed sneutrino VEV. In this subsection. 



we analyze the lepton violating interactions in (127). The question of baryon 



violation will be discussed in the next subsection. 

It is well-known |in]-|13] that the lepton number violating terms in (127) 
infiuence the baryon asymmetry at high temperature in the early universe. 
The requirement that the existing baryon asymmetry is not erased before 
the electroweak phase transition typically implies [6Tj that 



.10-6/ V 10 y ^ ' 

Parameter €3 for a given tan /3 can be explicitly evaluated for any B-L MSSM 



vacuum using (125). For example, consider the vacuum specified by point 
(P) in Figure 1. This has the values tan/3 = 18 and c^(0) = 1.0. RG running 
c^ down to the electroweak scale, we find that c^(tEw) = 0.855 and, hence, 
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that /i = 0.855 A^. The VEV of z/3 can be obtained using (107). For the 



parameters of this vacuum, (z/3) = 4.433 M.. Finally, unless otherwise stated 
we will take the third family neutrino Yukawa coupling to be X,^^ ~ 10"^°. 
This choice is motivated by the constraints on proton decay and will be 



discussed in the following subsection. Putting these values into (125) gives 
63 ~ 5.185 X 10^^° and, hence, 

(_£!,) (^). 0.33x10-3. (130) 

well below the necessary bound of unity. If we sample over all five vacua 
(P),(Q),(R),(S),(T) specified above, we find that 

0-688 X 10-^ < (^) (^) < 1.04 X 10-^ , (131) 



in each case below the bound in (129). We conclude that our B-L MSSM 



theory satisfies the conditions for baryon asymmetry. 

As discussed in [611[62], theories with lepton number violating interactions 



of the form in (127) naturally solve many fundamental cosmological problems 
if the gravitino is the lightest supersymmetric partner (LSP). The lifetime of 
the gravitino is then found to be [61] 



2 



Assuming that the lightest neutralino is the next-to-lightest superparticle 
(NLSP), one finds that 

— Mi^)l^)"^(.oa^)"'- (-, 

These results are relevant to the B-L MSSM theory discussed in this 
paper. First, it is possible to choose parameters so that the gravitino is, 
indeed, the LSP. Second, as can be seen from the spectra presented in the 
previous section at five different points, the lightest standard model sparticle 
is always the neutralino N^. As an example, let us compute the lifetimes of 
the gravitino and the lightest neutralino at the point (P) in Figure 1. From 
Table 2, we see that m^o = 147 GeV. Hence, adjusting the gravitino mass 
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to be, say, ■013/2 = 80 GeV, makes it the LSP while N^ is the NLSP. Using 



this value for ^3/2 and (130), it then follows from (132) that 



T3/2 



3.45 X 10^^^ 



(134) 



Noting that the age of the universe is typically estimated to be 13.7 billion 
years, that is, 4.32 x 10^''' seconds, we see that the gravitino lifetime greatly 
exceed this. Hence, the gravitino is the primary candidate for dark matter. 
On the other hand, using m^o = 100 GeV and (130), we find from (133) 
that 

Tr,LSP ^ 1.77 X 10-^s , (135) 

much to short-lived to form dark matter. Let us extend these results by 
evaluating the LSP and NLSP lifetimes at the five points (P),(Q),(R),(S),(T) 
specified above. Choosing m3/2 to be 20 GeV lighter than the corresponding 



N^ mass, we find using (131) that 



29, 



and 



1.65 X lO^^s < r3/2 < 2.47 x W'^s 



1.45 X lO-'^s < TNLSP < 5.52 x 10"^ . 



(136) 



(137) 



We conclude that for a gravitino LSP, our B-L MSSM theory has a long- 
lived gravitino consistent with it being dark matter, as well as an NLSP 
which decays very rapidly. 

5.2 Baryon Number Violation 

Recall that since U{1)b-l contains matter parity, the dangerous lepton and 
baryon number violating interactions in (7) are disallowed in the high energy 
superpotential. At much lower scales, the B-L violating VEV (u^) can po- 
tentially re-introduce these terms. As discussed above, however, this VEV 
induces from the dimension four superpotential only the lepton number vi- 



olating interactions in (127). The baryon number violating Uidjd^ terms are 
not regenerated. Therefore, to this order, baryon number is conserved and 
the proton is completely stable. However, the superpotential can contain 
B-L invariant higher dimensional terms proportional to u^Uidjdk- When the 
sneutrino develops a non-zero VEV, this generates effective dimension four 
operators of the form 

X'ljf^Uidjdk , (138) 
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d 




Figure 13: Effective operators generated by the dimension 4 interactions 
X^ji^LiQjdk and X'^j^Uidjdk- When the external fields are light families, these 
graphs induce nucleon decay. The solid lines represent fermions while the 
dashed line represents scalar propagators. 



where 






Mr 



(139) 



lijk = —likj are dimensionless parameters, and Mc is the compactification 



scale which we loosely identify with M„ in ( 16 ). For proton decay, the relevant 
operators are 

KikUididk (140) 

with k restricted to A; = 2, 3. 

Lepton number violating terms of the form X'-^^LiQ jdk can combine with 
the baryon number violating interactions X'l^fMidjdk to produce the effective 



operators in Figure [T3| Generically, these operators can induce proton decay 
via several channels. For the specific B-L MSSM theory in this paper, how- 



ever, it follows from (127) that the relevant lepton number violating terms 
are restricted to 

(^3>^d,kL3Qkdk (141) 

with k = 2,3. Since the t~^ and iJ'^-meson masses exceed that of the proton, 
in our specific theory the only potential decay channel is p — )■ K~^ + z/3. 



We find from (139), (140) and (141) that the product of the dimensionless 



couplings inducing this decay is 



XX— esAd 27112 



Mr. 



(142) 



As discussed in 
imental bound if 



5] , this channel will be suppressed below the exper- 
A'A" < C(10~25) _ (^43) 
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In estimating this bound, we have taken the mass of the intermediate squark 
in Figure 13 to be of (9(1 TeV), corresponding to its derived values in Section 
4. Under what conditions is (|143) satisfied? To be concrete, let us compute 



product (142) at the point (P) in Figure 1. As discussed above, here tan/3 
18, fi = 0.855A^ and (z/3) = 4.433A^. Leaving, for a moment, X^.^ arbitrary, 
one obtains e^ ~ S.lSSAjyg. Using this and (16), we find that for the p — >■ 
K^ + U3 channel 

A'A" = e3Ad,27ii2^ = 6.89 x 10-i^A,3 7ii2 • (144) 



Assuming, for simplicity, that 7112 is of C(l), it follows that bound (143) will 
be satisfied by taking 

K, < 10-1° . (145) 

We arrive at a similar conclusion for each of the remaining four points 
(Q),(R),(S) and (T). This explains our choice of the upper bound X^,^ ~ 10~^° 
in the previous subsection. Of course, choosing X^^ < 10"^° and/or 7112 < 1 
will suppresses proton decay even further below the experimental bound. 
Finally, although we relegate a detailed discussion of B-L MSSM neutrino 



masses to a future publication, it follows from (154) in Appendix A that the 



choice Ajyg ~ 10 will generate a third family left-handed neutrino mass of 
order 

m^3 ~ i^^^iM! ^ io-6ey . (146) 



m,7 



NO 



In evaluating rriN-i, we have taken the MSSM neutralino mass m^o to be of 



0(100 GeV) corresponding to the derived values in Section 4. The mass ( 146 ) 
is consistent with an "inverted" hierarchy of neutrino masses, where the third 
family neutrino is the lightest. We emphasize that our entire analysis has 
been for the specific B-L MSSM theory with diagonal Yukawa interactions 
and spontaneous B-L breaking exclusively by an expectation value of z/3. For 
generic off-diagonal couplings and B-L breaking involving the remaining two 
sneutrinos, the analysis of baryon asymmetry, proton decay and the neutrino 
mass hierarchy can be considerably generalized. 

Appendix A - Mass Diagonalization 

The B-L MSSM theory considered in this paper is a minimal extension of the 
standard supersymmetric model. Hence, portions of the analysis of the mass 
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spectrum is similar to that of the MSSM; see for example |73]. However, the 
addition of right-handed neutrino multiplets, the extra U{1)b-l gauge factor, 
the associated new soft SUSY breaking terms and the masses/field mixings 
induced by the (z/3) vacuum expectation value substantially complicate these 
calculations. In this Appendix, we briefly present the analysis of the mass 
spectrum used in this paper, emphasizing the differences from the MSSM. 

Neutralinos: 

We begin by discussing the mass eigenstates of the neutralinos. Restrict- 
ing ourselves, for the time-being, to the MSSM portion of our theory, the 
neutralino mass matrix will be of the form 



/ Ml -^ sxm\ 

92{H) g2{H) 

y/2 y/2 



M2 

9Y^il) 92(H) 

9y{H) g2{H) 



-fi 



(147) 



9Y\n! g2\nj _ PI I 



= ~ 

in the basis (Ay,Avt'0,if ,H). In this matrix, the fi entries arise from the 
Higgsino portion of the /i term, the M, entries are from the associated soft 
breaking terms and the off-diagonal entries are derived from the gauge super- 
field couplings in the Higgs kinetic energies. The physical neutralinos are the 
eigenstates of this matrix, which we label as N^ for i = 1, .., 4. The associated 
masses are denoted m^o. Note that when restricted to the MSSM portion of 

i 

our theory, the masses of both the left- and right-chiral neutrino fermions, 
ijjNi and ipij. for i = 1, 2, 3 respectively, arise from pure Dirac interactions and 
do not couple with the neutralino mass matrix. 

The full B-L MSSM theory discussed in this paper, however, is more 
complicated. There is a fifth neutral gaugino, that is, the superpartner \b-l 
of the Ab-l boson. Although this field does not mix with the above MSSM 
neutralinos, it does mix with the third family right-handed neutrino through 
the 1/3 kinetic energy. It follows that a complete description requires extend- 

ing the above basis to {ipNs, i^u3, ^b-l, Ay, Xw", H , H^). For this basis, the 



55 



mass matrix takes the form 



/ 




















V2gB-L{l^3) 

M4 












Ml 









M2 








9y{H) gYJH) 








gYJH) g2{H) 

/2 



9y{H) 
V2 



V2 
92{H) 



92{H) 





V2 
92{H) 



V2 



-/^ 



V2 

/ 

(148) 

We note that we have dropped all terms in this matrix proportional to the 
third family left-handed sneutrino expectation value. This has been calcu- 
lated in EDI ESI and found to be 



{N,) 






V"^^3 



For the vacua in this paper, we find 

(iV3)«A.3(z/3) 



(149) 



(150) 



and, hence, such terms can be safely ignored. 

In the limit of vanishingly small X^^, matrix (148) splits into the MSSM 



neutralino matrix in (147), the 2x2 mass matrix 



V2gB-L{i^3) 

V2gB-L{jy3) M4 



(151) 



mixing the right-handed neutrino tpu^ and the B-L gaugino Xb-l and a 



vanishing left-handed neutrino mass. Diagonalizing (151) assuming M4 <^ 



2\f2gB-Lk^zi ■, which is the case for the vacua described in this paper, gives 

M4 , ., ,. M4 , 



m^ 



MAe_.(l- 



l^plgB-L^y^) 



'"^A' , =Mas_^( 



1 + 



I^PlgB-Lkvz) 

(152) 

where we have used (67). Note that for M4 — )■ 0, that is, the supersymmetric 

limit in this sector. 



m 



V-. 



■mv 



Ma 



(153) 
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Hence, in addition to the B-L vector boson and gaugino masses becoming 
identical, as they must, there is also a predicted degeneracy between the right- 
handed neutrino mass eigenstate and the Ab-l mass. These degeneracies 
split as M4 is turned on. Now consider \^^ small, but non-zero. This gives 
small corrections to the B-L gaugino and right-handed neutrino masses. In 
addition, it generates a non- vanishing left-handed neutrino mass p5l[M] . For 
the vacua described in this paper, this is well-approximated by 

^ (KMl , (154) 



m 



NO 



where m^o is a typical MSSM neutralino mass. 

For the actual calculations in this paper, however, we do not use any of 



these approximations. Rather, we consider the complete mass matrix (148) 
and diagonalize it numerically. 

Charginos: 

Now consider the mass matrix for the charginos. Note that the non-vanishing 
sneutrinos VEVs in the B-L MSSM theory lead to off-diagonal chargino mix- 
ing terms which are not present in the MSSM. However, these are propor- 
tional to either (A^s) or X^j^iy^) and, thus, are very small, far below the order 
of approximation used in this paper. It follows that the chargino mass matrix 
is very nearly that of the MSSM given, for example, in [73]. In this paper, the 
chargino masses are calculated by diagonalizing their MSSM mass matrix. 

Squarks and Sleptons: 

Next, let us analyze the squark and slepton mass matrices. These have super- 
symmetric contributions from both the D- and F-terms in the scalar poten- 
tial energy. First consider the D-terms. In the MSSM, left-handed squarks 
and sleptons get well-known D-term masses proportional to the square of 
the Higgs VEVs. However, as discussed in [59], in the B-L MSSM each 
squark and slepton gets an additional-and significant-D-term contribution 
to its squared mass proportional to gs-ii^s) ■ ^'^^ example, these contribu- 



tions to the left-handed squark and slepton squared masses appear in (66). 
In evaluating scalar masses in this paper, we include these additional B-L 
contributions. Note that there will also be D-term contributions to some 
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scalar masses arising from the left-handed sneutrino VEV (A^s). However, as 
discussed above, this is very small and will be ignored in our calculations. 

Now consider the F-term contributions. Recall that we have assumed in 
^ that the squark/slepton Yukawa couplings are diagonal and from ( |145 ) 



that A;^3 must be very small. It follows that the non-zero (z/3) and (N^) 
VEV's in the B-L MSSM do not induce significant contributions and, hence, 
the F-term scalar masses can be well-approximated using the pure MSSM. 
Note that since the Yukawa couplings have been chosen to be diagonal, the 
only possible mixing is between the left and right-handed states of each scalar 
field. For the "up" scalars 0^, this mixing arises from 

||^r = |/i^ + A<^C,L</'«,i?r (155) 

when the Higgs obtains a VEV. Similarly for the "down" type scalars. Since 
the Yukawa couplings for the third family are considerably larger than the 
others, it is reasonable to drop the first and second family contributions to 
the F-term scalar masses. 

Finally, note that there are also non-supersymmetric contributions from 
the quadratic and cubic soft terms. The quadratic masses are diagonal and 
straightforward. A contribution to mixing terms arises from the soft cubic 



couplings. Recall from (26) that the A-parameters are assumed to be pro- 
portional to the corresponding Yukawa couplings. Hence, we drop all such 
mixings except for the third family squarks and sleptons. 

Putting everything together, in the basis (0l, (Pr) the squark and down- 
slepton mass squared matrices m^, m^, and m^ are given by 



^2_ ( rnl^ + mi 


m\ + m\ ) 


(156) 




{H)Al-f,X,{H)\ 
ml+ml J 


(157) 


2 / m|g + ml 


{H)A;-fiXAH)\ 
ml +ml J- 


(158) 



These can be diagonalized by a unitary matrix. For example, in the case 
4> = t^ this unitary matrix can be written as 




c, -s*A (tj 



hi C-i j \ t/O 
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(159) 



To obtain the masses m?(i)m?(2) for the eigenstates t^ ,t^ , matrix mjin^ is 

diagonahzed and one takes the square root of the eigenvalues. However, if 
any of these are degenerate, the Takachi diagonahzation process ^6] must 
be used. The sneutrino mass squares were analyzed in detail in [59]. Hence, 
we do not discuss the up-slepton mass matrix in this paper. 

Higgs Bosons: 

Next, consider the masses of the five Higgs bosons. We note that scalar po- 
tential is perturbed from the MSSM potential by the VEVs of both the right 
handed and left handed neutrinos. However, similar to previous discussions, 
we note that these perturbations are small and choose to approximate it with 
the MSSM scalar potential. 

As is well known, in the MSSM the tree level mass of the lightest neu- 
tral Higgs is bounded above by the mass of the Z boson [73]. However, a 
neutral Higgs boson with mass below Mz in the MSSM has been ruled out 
experimentally. Hence, one must explore the radiative corrections to Higgs 
masses. It turns out these corrections to the lightest neutral Higgs mass are 
quite sizable and allow the MSSM to remain a viable theory, albeit with a 
restricted range of parameters. 

We note that there exists an extensive literature on the radiative correc- 
tions to the MSSM Higgs masses, with leading corrections calculated up to 
the three-loop level [HZ]. For the purposes of this paper, we consider the 
lightest neutral Higgs mass to the first few leading terms in the one-loop 
correction, as given in |88j . That is, our approximation for the lightest Higgs 
mass is 

'^lo = -{m\o + Ml + Ut) 



m 



^,Ml cos2 2/3 + ^ (m^o - M^ 



(160) 
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where 



Q , m-{i)m-(2) {^%) ~ ^%)) ^%) 

to 



■-3 



_i_ „4 4 ^^3 '^a I 1 — — /' ^ 



"tl t2 ,m^ \ 9^m2 



2 

<1) ,''^,-(2) 



m^ V 2 m_(2) — ""^ 



3 '•a 



rti) ■■-,-12) N >, 

^)ln H^)|- (161) 



The one-loop corrections to the remaining neutral and charged Higgs masses 
are also presented in [88j. It is useful, however, to note two things; first, 
the one-loop corrections to the A^ mass are sub-leading when inserted into 
(160) and, second, the radiative corrections to the A^ fl'^^^ and H^ masses do 



not substantially change conclusions about their lower bounds derived from 
their tree-level expressions. That is, the masses for the remaining neutral 
and charged Higgs can be well-approximated by their tree level expressions. 
These are given by [HE] 



m^o = 


= 2/ip + m^ + m^ 


(162) 


m^o = 


- i«.+Mi).;«»;^'i'^ 


-m^oM|cos2(2/3)) (163) 


m^± = 


= m^o + M^ . 


(164) 



We have followed the standard notation of labeling the Higgs mass eigen- 
states as /i°(iJ°) for the lightest (heaviest) neutral Higgs, J^ for the CP odd 
neutral Higgs and H^ for the respective ±-charged Higgs. It is also useful 
to recall that 

Ml = \glmf ^ {Ef) . (165) 



The definition for Mz was given in (28) and (59). Finally, note from (|160) 



and (163) that 



muo > m/jO (166) 

at tree-level and using [88] that this inequality persists when the one-loop 
corrections are included. 
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Appendix B - Experimental Bounds on Higgs 
Masses 

In this Appendix, we review the experimental bounds on the Higgs masses. 
We point out that these bounds are highly model dependent and, even within 
the MSSM, depend on assumptions imposed on the soft SUSY breaking pa- 
rameters, the structure of mass matrices and so on. To illustrate the non- 
trivial nature of these bounds, we briefly review one scenario of the MSSM 
and the experimental results for the Higgs masses. Since Higgs fields carry no 
B-L charge, to the one-loop level the following discussion is directly relevant 
to our U{1)b-l MSSM theory. 

We will discuss what is commonly referred to as the myfi — max scenario, 
in which one assumes CP conservation and that the off-diagonal terms in the 
stop mixing matrix are large. For this scenario, the Higgs mass bounds are 



reported in ^75l l89j and reproduced here in Figure 14 We begin our analysis 



with plot (c) of Figure 14, which presents the experimentally excluded zone 



in the m^o-tan/3 plane. Recall from (31 ) that to guarantee that the Yukawa 



parameters are perturbative, we have restricted the range of tan/3 to be 

4<tan/3<50. (167) 

Note from (c) that the lower bound on uiaq can be as small as ~ 93 GeV, but 
only for large values of tan/3. Furthermore, for 93 GeV < niAo < 120 GeV 
the lower bound on tan/3 must exceed 4 in a complicated way. Therefore, to 
simplify the calculations in this paper we will always take 

ruAO > 120 GeV . (168) 

With this assumption, there is no restriction on the the value of tan/3 beyond 



those of (167). 



We see from plots (a) and (b) respectively that the lower bounds on 



the lightest Higgs mass m^o, and, hence, using (166), on ruffo, are highly 



dependent on the values of m^o and tan/3. First, consider plot (a). For 

rriAO > 120 GeV =^ rriho > 114.4 GeV , (169) 

the oft-quoted lower bound on the lightest neutral Higgs mass. Note that 
for 93 GeV < m^o < 120 GeV there is a small, almost linear, allowed region 
where 94 GeV < rriho < 114.4 GeV. Theoretically, finding parameters in this 
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r igUre 14; Exclusions at the 95% confidence level (medium-grey or light-green) and at the 99.7% 
confidence level (dark-grey or dark-green) for the CP conserving, nif^o — max scenario and for nit = 174.3 
GeV/c^. The figure shows the theoretically inaccessible domains (light-grey or yellow) and the regions 
excluded in this search, in four projections of the MSSM parameters: (a): (m^jO, »?i^o); (b): (m^o, tan/3); 
(c): (m^o, tan/3); (d): {mjf±, tan/3). The dashed lines indicate the boundaries of the region which are 
expected to be excluded at 95% confidence level based on Monte Carlo simulations with no signal. In the 
(m^o, tan/3) projection (plot (b)), the upper boundary of the parameter space is indicated for four values 
of the top quark mass; from left to right: mt =169.3, 174.3, 179.3, and 183.0 GeV/c^. 
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region would requiring an unnaturally high degree of fine-tuning. This is 



another justification for choosing the bound in (168). Second, let us analyze 



the bounds in plot (b). It follows from (167) and (169) that plot (b) puts no 
further restrictions on the parameters. 

Now consider the remaining plot (d) in Figure 14. This presents the 



experimental exclusion zone for raH± as a function of tan/3. Note from (164) 



that when m^o exceeds its lower bound in (168), then mui must be larger 
than ~ 144 GeV. That is, 

m^o > 120Ge1/ =^ mH± > lUGeV . (170) 

With the exception of the small region 144 GeV < mH± ^ 160 GeV, where. 



roughly, tan/3 must exceed 6, the bounds (167) and (170) guarantee that plot 



(d) puts no further constraints on the parameters. 

To conclude: in the numerical analysis in this paper we use the mf^o —max 
scenario to put lower bounds on all five Higgs masses. This is accomplished 



by first calculating m^o using (162), accepting the result if m^o > 120 GeV 
and rejecting it if rriy^o has a smaller value. Assuming tti^o satisfies this 



inequality, we then calculate rriho from (160) and (161). Note that to evaluate 



cut in (161), one must specify a value for the top quark mass. In order to 



compare our results to Figure 14, we take 

mt = 174.3 GeV . 



(171) 



If 



The result is accepted if rriho > 114.4 GeV and is rejected otherwise. 
rriho exceeds this bound we continue by computing mH± using (164). If it 



lies in the range 144 GeV < mH± < 160 GeV, then tan/3 is restricted to 
exceed 6. If it does we accept the result and reject it if tan/3 is below this 
value. For r7i//± > 160 GeV, there are no further constraints and the result is 
accepted. In this way, for a given point in the Cq{0)-Ci,^{0) projection we can 
map out the allowed Higgs scalar region in the c^(0)-tan^ plane. Outside of 
this region, at least one of the five Higgs masses is less than the experimental 
lower bound. 

For a more thorough discussion of this and many other scenarios in the 
MSSM, see 



Appendix C - Kinetic Mixing 



As is well-known, see for example [BSlEn], in theories with at least two f/(l) 
gauge factors mixing can occur between the kinetic terms of the associated 
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vector supermultiplets. More specifically, in the B-L MSSM theory discussed 
in this paper there can be mixing between the U{1)y and U{1)b-l gauge 
kinetic terms. To analyze this, consider the hypercharge and B-L generators 
Y/2 and Yb-l defined in Section 2. Then 

Tr{jf = 11, TrYBV = 16, Tt^Yb-l = 8 . (172) 

The fact that Tr^Y^-h 7^ implies that kinetic mixing does indeed occur, 
introducing a new dimensionless coupling for the mixed term which we denote 
by gM- Following [68], we construct the 2x2 matrix gai3 with components 
9ii = 9y, 922 = 9B-L, 9i2 = 921 = 9m and define the new generators 

2 
Y^ = J2^'0Pc., a = 1,2 (173) 

/3=1 

where Y'^ = Y/2 and Yg = Yb-l- It then follows that at the one-loop level 
the RGEs for the gauge parameters generalize to 

^ = Y^/3a/^, (3^p=9a,Tr%Y^. (174) 

The expressions for /3q^ are cubic polynomials in the gauge couplings and 
easily worked out using Y/2 and Yb-l- If one assumes that qm ^ 9y, 9b-l, 
then 

/3n = Py = n9Y, /322 = Pb-l = lQ9l^L, /3i2 = Pm = 8glgB-L ■ (175) 

Note that the first two /3-functions are exactly the ones used to scale gY and 
gB-L in this paper. However, it is important to be cognizant of the existence 
and scaling of gM and its physical implications. 

Using the scaling solutions for gY and gB-L, one can solve the RGE for 
gM- The result is 

/ N 4^(0) f , 1 / 33^(0)2|tl 
guKt) = ^m(O) y=- I arctan(- yil + — ) - arctan( 

(176) 




with g{0) given in (15). Note that 5'm(^) decreases monotonically as one scales 
down from the unification scale t = to the electroweak scale tEw — —33.3. 
We find that 

gMitEw) = 9MiO) - .U89 . (177) 
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The initial value gM^O) is a model and threshold dependent quantity. As with 
the soft SUSY breaking parameters, its calculation requires specifying an 
explicit ultraviolet vacuum-which is beyond the scope of this paper. Hence, 
as with the soft parameters, we treat it as an arbitrary parameter. To make 
the following discussion explicit, we find it convenient to choose gM^^) = 
.1489/2. It follows that gM has the positive value .0745 at the unification 
scale, passes through zero at lower energy and becomes gnitEw) = —.0745 
at the electroweak scale. Over this entire range the assumption that gM ^ 
gY,gB-L remains valid. 



What is the physical consequence of ghiitEw) 7^ 0? It follows from (173) 
that, in addition to the hypercharge generator gri^), the mixing induces a 
new generator 

g^{a^ + bYB^L), (178) 

where 

gxa = gM, gxb = gM + gB-L ■ (179) 

It was shown in [64J that the parameter controlling the physical Z-Z' mixing 
is given to leading order by 






2i = --y^i-^^y I v^i ^--v^i I . (180) 



Using the data in this paper and the choice (7m(0) = .0745, this becomes 



M. ^ ^ 



2^ ~ .5292 ,^ . (181) 

The hierarchy Mz/Mab_l is sufficiently large everywhere in the allowed re- 
gions discussed in Section 4 that 

2^ < 0(10-=^) , (182) 

which is phenomenologically acceptable |Mj. For example, at the points (P) 



and (S), where Ma^^^ = 1252 GeV and Ma^.^ = 3005 GeV, (|l8l|) becomes 
2^ = 2.796 X 10"^ and 2^ = 4.853 x lO""^ respectively. Finally, we note that 
taking larger f^M (0) will reduce the mixing at the electroweak scale, with 2^ 
vanishing for (^a/ (0) = .1489. On the other hand, for smaller values of 5'a/(0), 



such as f7A/(0) = 0, the Z-Z' mixing continues to satisfy (182) over most of 
each allowed region in Section 4. 
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We conclude that there is a reasonable range of (^Af (0) for which 1) the 
scaling equations for the gauge parameters are well-approximated by drop- 
ping the qm contributions and 2) for which the induced Z-Z' mixing is natu- 
rally small and phenomenologically acceptable. The analysis in this paper is 
meant to establish the existence of 5-L/electroweak symmetry breaking with 
a reasonable spectrum over a subspace of a large initial parameter space, and 
not to give precision values of all quantities. It is consistent with this goal, 
and greatly simplifies the calculations, to ignore the U{1)y-U{1)b-l kinetic 
mixing. 
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